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ABSTRACT  B 

This  study  investigated  the  problem  of  synthesizing  a minimum  cost 
transportation  system  plan  to  service  forecast  shipment  requirements  among 
a set  of  points  in  such  a manner  as  to  satisfy  aggregate  ship-time  performance 
levels.  Both  commercial  and  dedicated  modes  may  be  used  in  the  transporta- 
tion plan;  but  the  latter  must  be  designed  in  detail  including  specification  of  ve- 
hicles to  be  used,  the  route  each  is  to  service,  and  arrival/departure  time 
schedules,,  The  problem,  which  is  important  to  operations  of  the  Air  Force 
Logistics  Command,  requires  prescription  of  a transport  network  capable  of 
providing  acceptable  levels  of  service  to  the  shipments  which  it  accommodates. 

Research  objectives  included:  refinement  of  a modeling  approach  ini- 
tiated during  the  1976  USAF/ASEE  Summer  Faculty  Program,  development  of 
solution  approaches,  and  assessment  of  computation  time  necessary  to  solve 
problems  of  realistic  size.  Study  results  satisfy  each  objective. 

A large-scale  mixed,  01  integer,  linear  programming  model  of  the  plan- 
ning problem  is  developed  and  simplified  for  solution  by  applying  Benders'  de- 
composition to  yield  two  more  simple,  interacting  subproblems.  One  of  these, 
a linear  program  which  assures  ship-time  performance,  is  amenable  to  efficient 
solution  techniques  (Generalized  Upper  Bounding  and  column  generation)  for 
which  specialized  algorithms  are  presented. 

Several  formulations  of  the  other  subproblem,  which  defines thododieatod 


lode  network,  were  developed.  The  first  remiires  enumeration  of  a set  of 
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feasible  vehicle  tours  from  which  the  dedicated  mode  may  be  designed.  A 
special  purpose,  implicit  enumeration  algorithm  applicable  to  this  model  is 
described. 

Three  additional  formulations,  each  of  which  constructs  vehicle  tours 
directly,  were  investigated.  One  model,  a linear  program  amenable  to  large- 
scale  programming  techniques,  constructs  tours  via  a column  generation  pro- 
cedure and  may  be  used  in  a branch-and-bound  solution  algorithm.  An  alternate 
model  casts  the  problem  as  a modified  material  flow-circulation  problem.  Char- 
acteristics of  this  model  which  may  be  exploited  to  devise  an  efficient,  implicit 
enumeration  algorithm  are  described.  Finally,  approaches  based  on  constraint 
aggregation  and  Generalized  Lagrange  Multipliers  were  investigated  in  the  at- 
tempt to  reduce  the  tour  construction  problem  to  a straightforward  flow- 
circulation  problem. 

In  all  cases,  the  tactical,  vehicle  scheduling  problem  is  treated  subse- 
quently using  an  available  linear  programming  approach.  Collectively,  the  re- 
sults appear  to  offer  the  capability  to  solve  transportation  system  planning  prob- 
lems of  realistic  size. 
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CHAPTER  I 


INTRODUCTION 

The  Air  Force  operates  a complex  logistics  system  which  is  composed 
of  transportation,  inventory,  and  repair  subsystems.  These  components  are 
related  by  the  function  each  performs  and  through  the  management  policies 
which  guide  system  operations.  The  transportation  subsystem  is  of  particular 
importance,  since  it  has  a significant  influence  over  the  capital  investment  re- 
quired for  spare  parts  and  since  it  provides  a critical  service  necessary  for 
maintaining  acceptable  numbers  of  weapons  systems  in  operational  condition. 
This  study  is  directed  to  the  design  of  the  transportation  system  operated  by 
the  Air  Force  Logistics  Command  (AFLC).  However,  the  mathematical  struc- 
ture of  the  problem  is  common  to  a variety  of  systems  (including  multi-modal 
transportation  and  other  collection-distribution  networks  such  as  that  operated 
by  the  postal  service)  and  may  be  specialized  for  application  in  other  contexts. 

In  order  to  optimize  the  entire  logistics  system,  an  optimal  balance 
must  be  achieved  amongst  its  operating  components:  transportation,  inventory, 
and  repair.  If  resupply  time  be  reduced,  fewer  investment  items  need  be  pur- 
chased to  provide  equivalent  logistics  support  to  accomplish  the  prescribed  Air 
Force  mission.  Since  many  investment  items  managed  by  the  Air  Force  are 
very  costly,  it  is  important  to  be  able  to  prescribe  an  optimal  trade-off  between 
inventory  investment  and  transportation  expense.  Dr.  A.  Kovacs,  a colleague 
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on  the  1 97(5  USAF/ASEE  Summer  Faculty  Research  Program,  has  examined 


the  logistics  system  problem.  He  developed  several  models  which  describe 
essential  features  of  the  logistics  system  and  showed  how  they  might  be  used  in 
a prescriptive  model  to  define  the  optimal  trade-off  between  transportation  ex- 
pense and  inventory  investment.  In  order  to  accomplish  this  overall  objective, 
Dr.  Kovacs'  logistics  system  model  requires  the  functional  relationship  of 
minimum  transportation  cost  to  service  level.  The  Transportation  System 
Planning  Model,  which  is  the  topic  of  this  research,  can  define  that  necessary 
functional  relationship.  Results  of  this  study  may  therefore  be  used  to  optimize 
the  transportation  system  plan  as  well  as  to  provide  guidance  in  making  impor- 
tant decisions  which  integrate  operations  of  the  entire  logistics  system. 

Problem  Statement 

The  essential  features  of  the  transportation  system  planning  problem 
are  stated  concisely  below  so  that  all  components  of  the  problem  addressed  are 
identified. 

Given  expected  daily  shipping  requirements  amongst  all  points  in  the 
logistics  system  categorized  by  the  shipping  times  required  and, 
perhaps,  by  characteristics  which  might  be  used  in  vehicle/mode 
assignment  decisions  (i.  e. , weight,  volume,  special  handling 
needs). 

Design  an  optimal,  minimum  total  annual  cost  transportation  system 
plan  which  defines  specific  routes  by  which  shipments  are  to  be 
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made,  assures  required  si  pping  time  for  each  material  category, 
and  determines  the  appropriate  balance  in  the  use  of  commercial 
modes  and  dedicated  modes,  which  should  be  specified  if  they  prove 
to  be  cost  effective. 

To  provide  the  necessary  flexibility,  the  dedicated  modal  network  could 
(potentially)  use  a number  of  different  types  of  aircraft  and  trucks. 
The  optimal  system  design  must  specify  dedicated  modes  in  de- 
tail including  the  best  types  and  numbers  of  vehicles  to  be  used 
and  designation  of  the  route  to  be  served  by  each  vehicle.  Char- 
acteristics of  each  vehicle  such  as  weight/volume  capacity,  speed, 
endurance,  and  maintenance  positioning  requirements  need  be 
considered.  The  arrival/departure  schedules  of  each  vehicle  at 
each  point  it  services  must  be  specified  since  this  schedule  is 
important  to  coordinating  dock  operations,  to  vehicle  crew  needs, 
and  to  service  levels  afforted  material  flows. 

The  system  problem  encompasses  three  important  aspects  of  design:  vehicle 
route  synthesis,  material  flowtime  control,  and  determination  of  vehicle  sched- 
ules. The  optimization  model  should  permit  detailed  sensitivity  analyses  of 
important  factors  such  as  shipping  time  requirements  and  transportation  prior- 
ity scheme.  The  model  should  also  facilitate  evaluation  of  management  policies 
such  as  the  use  of  a single  point  to  slock  second-echelon  inventories.  Develop- 
ment of  these  post-optimality  analysis  procedures  is  not  a part  of  this  study. 
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but  the  need  for  such  detailed  evaluation  was  identified  and  used  to  guide  de- 


velopment of  the  basic  modeling  capability. 

A brief  description  of  the  operating  procedures  used  in  the  current  sys- 
tem is  given  in  Appendix  A.  The  description  indicates  the  practical  importance 
of  each  facet  of  the  problem  statement  and  its  relationship  to  AFLC  needs. 

Summary  of  Approach  and  Report  Outline 
A model  of  the  transportation  system  planning  problem  was  developed 
during  the  1976  USAF/ASEE  Summer  Faculty  Program  and  presented  in  the 
proposal  submitted  to  AFOSR  which  lead  to  this  study.  The  original  model, 
which  is  given  in  Appendix  B of  this  report,  incorporated  the  essential  features 
of  the  transportation  system  planning  problem  but  offers  little  opportunity  to 
develop  efficient  solution  procedures  because  of  the  large  number  of  integer 
variables  and  constraints  in  the  formulation.  A significant  improvement  in 
formulation  was  developed  and  submitted  to  AFOSR  while  the  research  proposal 
was  in  review.  The  second  model,  which  is  listed  in  Appendix  C,  employed  a 
modeling  approach  which  allowed  a significant  reduction  in  the  number  of  vari- 
ables and  constraints  required  to  express  the  problem  and  offered  a better  op- 
portunity to  develop  successful  solution  approaches. 

Research  during  this  study  has  provided  further,  significant  improve- 
ments in  model  formulation,  resulting  in  the  forms  presented  in  this  report. 
Experience  gained  in  developing  the  first  two  formulations  Iod  to  the  more  re- 
cent version  which  is  much  more  simple  and  which  is  amenable  to  efficient  so- 
lution procedures  which  were  developed  during  this  study. 

L 
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Solution  strategy  allows  the  vehicle  scheduling  component  of  the  prob- 
lem to  be  treated  separately  according  to  a method  developed  by  a study  con- 
ducted at  the  Air  Force  Institute  of  Technology  (A FIT).  The  material  flow  and 
vehicle  routing  components  are  treated  jointly  in  each  of  several  solution  pro- 
cedures. 

Given  a network  of  tours  which  are  traversed  by  vehicles  in  dedicated 
modes  and  knowledge  of  commercial  modes  available,  the  planning  problem  re- 
duces to  a multi-commodity  material  flow  problem  subject  to  the  important 
ship  time  constraints.  A model  representing  this  flow  problem  is  developed  in 
Chapter  III  and  efficient  solution  procedures  are  described.  The  model  is  a 
large-scale  linear  program  which  is  amenable  to  Generalized  Upper  Bounding 
and  column  generation  techniques.  This  model  might  be  used  to  evaluate  a pro- 
posed, dedicated  mode  network,  or,  perhaps,  to  design  such  a network  heuris- 
tically.  Most  importantly,  the  model  is  used  as  a basic  component  of  more 
comprehensive  models  which  also  incorporate  the  vehicle  routing  component. 

Vehicles  made  available  for  use  by  the  dedicated  mode  network  are 
assumed  to  be  defined  by  applicable  operating  characteristics  such  as  carrying 
capacity,  tour  length  limitations,  and  tour  origin/termination  points  (which  are 
dictated  by  vehicle  maintenance-facility  location).  One  approach,  which  was 
developed  to  provide  an  optimal  transportation  system  plan  relies  upon  the 
assumption  that, at  least  in  some  applications,  it  is  efficient  to  define  a set  of 
feasible  vehicle  tours  from  which  the  optimal  dedicated  mode  network  may  be 
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designed. 


This  approach,  which  is  described  in  Chapter  IV,  may  use  either  a 
limited  set  of  tours  which  is  likely  to  include  the  optimal  solution,  or  the  set 
of  all  feasible  tours  which  may  be  enumerated  rapidly  for  problems  of  limited 
size.  A large-scale  mixed,  01  integer,  linear  programming  model  is  shown  to 
include  both  the  material  flow  and  vehicle  tour  components.  Efficient  solution 
procedures  are  developed  for  this  model  applying  Benders'  decomposition  to  de- 
rive two  interacting  subproblems  which  invoke  existing  theory  to  determine  op- 
timal transportation  system  plans.  One  of  the  subproblems  is  (essentially)  the 
same  as  the  material  flow  model  discussed  in  Chapter  III  and  the  other  is  ame- 
nable to  efficient  solution  by  a special  purpose,  implicit  enumeration  algorithm 
designed  specifically  for  that  purpose. 

A third  — and  more  highly  sophisticated  approach  — was  researched 
and  is  described  in  Chapter  V.  This  approach  relies  upon  the  assumption  that, 
in  some  applications,  it  may  be  best  to  employ  a solution  procedure  which  is 
capable  of  constructing  vehicle  tours,  rather  than  merely  selecting  tours  from 
among  a set  which  is  provided  for  use.  Three  strategies  are  described  to  im- 
plement this  approach  by  which  vehicle  tours  are  constructed  as  an  integral  com- 
ponent of  the  solution  procedure. 

The  first  strategy  leads  to  a large-scale,  mixed  integer  linear  program- 
ming model.  The  program  is  simplified  by  application  of  Benders'  decomposition 
yielding  the  material  flow  (linear)  program  and  an  interacting  integer  problem 
which  completely  defines  the  vehicle  network.  This  binary  problem  bears  a 
close  relationship  to  well  known  network  problems  of  the  flow  circulation  type. 
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Characteristics  which  allow  efficient  solution  and  techniques  to  exploit  these 
characteristics  in  an  implicit  enumeration  algorithm  are  discussed. 

The  second  strategy  gives  rise  to  a model  formulation  which  is  amenable 
to  efficient  branch-and  bound  procedures.  This  strategy  incorporates  linear 
programming  techniques  such  as  Generalized  Upper  Bounding  and  column  genera- 
tion using  a shortest  path  subproblem. 

The  third  strategy  examines  the  opportunity  to  construct  vehicle  tours 
using  an  efficient  flow  circulation  algorithm.  Techniques  for  aggregating  con- 
straints into  the  objective  to  yield  a pure  flow-circulation  problem  are  evaluated. 
Generalized  Lagrange  multiplier  techniques  are  also  considered  as  a means  of 
effecting  this  more  simple  structure. 

A literature  review  which  describes  prior  work  in  each  of  the  three  areas 
related  in  the  transportation  system  planning  problem  appears  in  Chapter  II.  In 
particular,  the  vehicle  scheduling  approach  developed  in  the  A FIT  study  is 
described  in  that  chapter. 

A final  chapter  records  conclusions  from  this  study  and  outlines  recom- 
mendations for  continued  research  on  the  transportation  system  planning  prob- 
lem. A large  number  of  variables  interact  in  the  planning  problem  and  the  no- 
tation required  to  represent  the  problem  is  lengthy.  A complete  listing  of  all 
notation  used  in  the  body  of  this  report  is  therefore  given  in  Appendi>  D for 
convenience. 
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CHAPTER  II 
LITERATURE  REVIEW 

A rather  thorough  search  of  the  literature  has  revealed  that  there  is 
little  history  on  the  stated  problem  which  requires  simultaneous  optimization 
of  a network  design  and  control  of  multicommodity  material  flowtimes.  There 
are,  however,  several  recent  studies  which  are  of  interest  because  of  the  ap- 
proaches they  take  to  solving  related  problems.  There  is  also  a well  developed 
literature  on  each  of  the  three  main  aspects  of  the  planning  problem — vehicle 
network  design,  multicommodity  material  flows  through  networks,  and  vehicle 
scheduling.  Each  of  these  areas  is  reviewed  briefly  in  this  section. 


Related  Problems 

Two  master's  theses  addressed  to  the  design  of  the  Logair  route  struc- 
ture have  recently  been  completed  at  the  Air  Force  Institute  of  Technology 
(A FIT).  Building  on  the  experience  gained  by  the  Palmatier  and  Prescott  (1975) 
study,  McPherson  and  O'Hara  (1976)  formulate  a limited  version  of  the  network 
flow  design  problem  which  minimizes  a total  cost  function 

min  2 (dC,.  + f)  Z.. 
ij  iJ  ' U 

subject  to  range  limitation  on  sortee  length,  material  conservation  constraints 

2 X...  + R..  = 2 X...  , V .1, 

ij  ijk  jk  £ jlk  Ifk 

i + j.  i fk 


vehicle  capacity  restrictions 
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S X...  £ d..  Z..  ; V.  1 1 
k ijk  ij  1J  ifk 


and  vehicle  conservation  constraints 


in  which 


2 Z..  - 2 Z.. 
i i J1 


X...  s 0 ; Z..  = 0, 1, . . . ,n 
ijk  i J 


d = cost  per  mile 

C.  - distance  from  point  i to  j 
ij 

f aircraft  landing  fee 

Z - number  of  times  arc  ij  is  flown  = 0, 1,  2, ....  n 
X..^  = tons  shipped  from  i to  j bound  for  destination  k 
d„  = aircraft  capacity  on  arc  ij 
R.,  = total  requirements  originating  at  j for  k 

This  model  is  a large  scale,  mixed  integer  linear  program  (MILP)  which  the 
authors  solved  using  a standard,  although  sophisticated,  algorithm  available  on 
the  AFLC  computer  system.  Three  serious  limitations  arise  in  this  formulation. 
First,  optimal  tours  are  not  specified  for  vehicles.  Secondly,  transportation 
time  for  a particular  shipment  cannot  be  determined  because  the  model  considers 
only  total  cargo  shipped  out  of  each  point  rather  than  accounting  for  the  route  of 
each  shipment.  This  omission  will  not  allow  a useful  measure  of  service  level 
to  be  developed.  Finally,  it  took  three  hours  of  computation  time  (on  a Honey- 
well fif>0  computer)  to  solve  a twelve  point  problem,  so  applicability  of  this 
approach  is  limited  to  small  network  problems. 
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Three  other  studies  are  of  interest  because  they  successfully  apply  large 
scale  programming  techniques  to  related  network  design  problems.  Rao  and 
Zoints  (1968)  studied  the  problem  of  minimizing  the  cost  of  transporting  com- 
modities amongst  various  ports  by  vessels  which  are  routed  according  to  model 
prescription.  They  employ  a column  generation  scheme  for  which  subproblems 
are  solved  by  the  efficient  out-of-kilter  algorithm.  Computationally,  the  approach 
is  efficient,  but  integer  solutions  do  not  result  and  connected  vessel  routes  are 
not  specified,  so  application  of  the  model  is  limited. 

Richardson  (1976)  formulated  an  airline  routing  design  problem  as  a 
MILP  problem  and  solved  it  by  application  of  Benders'  decomposition  method. 

A standard  linear  programming  "subproblem"  and  a 0,  1 integer  "master"  prob- 
lem resulted  and  the  latter  was  solved  by  a special  purpose  algorithm.  Compu- 
tation times  for  problems  of  realistic  size  (26  points)  were  modest  (less  than 
two  minutes  on  a DEC-10/1055  computer). 

Geoffrion  and  Graves  (1974)  studied  a multicommodity  distribution  prob- 
lem involving  14  manufacturing  plants,  45  possible  distribution  center  sites, 

121  customer  zones  and  17  commodity  classes.  They  solved  the  large  scale 
MILP  problem  by  a Benders'  decomposition  formulation  which  efficiently  solved 
subproblems  of  the  transportation  type  and  a master  0, 1 problem  using  a speci- 
ally designed,  hybrid  branch-and-  bound/cutting  plane  approach.  Computation 
times  less  than  one  minute  using  an  IBM  360/91  computer  were  reported.  The 
authors  also  present  formidable  arguments  supporting  their  optimization  ap- 
proach since  it  allows  comprehensive  sensitivity  studies  and  detailed  analysis 


of  managerial  alternatives. 
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While  the  last  three  studies  were  not  addressed  specifically  to  the  trans- 
portation system  problem,  they  do  indicate  the  computational  success  of  certain 
large  scale  programming  techniques  in  solving  MILP  models  of  related  form. 
Other  examples  could  be  cited,  but  these  should  provide  sufficient  background 
for  later  discussion. 


Vehicle  Network  Design 

The  classical  approach  to  designing  a network  through  which  vehicles 
transport  goods  is  through  the  multiple  traveling  salesman  formulation.  This 
problem  requires  each  of  n points  to  be  serviced  once  by  only  one  of  the  m vehi- 
cles provided  in  the  system;  vehicles  must  collectively  execute  their  itinerancy 
at  minimum  total  cost.  This  problem  may  be  expressed  mathematically  (see 
Taha  (1971,  Svestka  (1973),  and  Miller,  et  al.  (i960)  ) as: 


min  £ d..  x.. 
ij  i]  iJ 

[11 

subject  to: 

S X =1 

j=l,  2 n 

[21 

i lJ 

£ X..  = 1 

j lJ 

i=l,  2, . . . , n 

[31 

Y.  - Yj  + (n-m) 

Xjj  s n-m-1 

[41 

in  which  d^  is  the  cost  of  traversing  the  arc  from  i to  j and  X^  is  1 if  a vehicle 
travels  over  the  arc  between  points  i and  j and  "0"  otherwise.  The  last  con- 
straint requires  each  vehicle  tour  to  be  a single,  connected  cycle  in  the  network. 
Different  formulations  to  invoke  the  connected- cycle  requirement  have  given  rise 
to  a number  of  mathematical  models  which  define  the  basic  traveling  salesman 
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problem.  Golden,  Magnanti,  and  Nguyen  (1977)  provide  a rather  complete  tax- 
onomy of  integer  formulations  of  this  type. 

Balinski  and  Quandt  (1964)  suggested  an  alternate  type  of  program  which 
is  essentially  a set  partitioning  problem: 


min 


J 


2 Cj  X. 

1=1  J J 


[51 


s.  t. 


£ a(j  Xj 


J=1 


= l; 


i=l.2. 


■ m 


V°-i' 


j~l.  2, . . . , n 


[61 


[71 


in  which 


n is  the  total  number  of  tours  satisfying  vehicle  capacity  constraints. 
a„  = l if  delivery  route  j visits  point  i,  0 otherwise 
Cj  = cost  of  delivery  route  j 
Xj  = 1 if  delivery  route  j is  used,  0 otherwise. 

The  authors  developed  a column  generation  technique  to  solve  problem  [5] 
through  [71.  The  approach  has  apparently  not  been  further  refined  by  more  re- 
cent research. 

The  (multiple-)  vehicle  routing  problem  has  received  a great  deal  of 
attention  and  a number  of  solution  approaches  have  been  developed.  Three  of 
these  general  approaches  are  of  interest  with  respect  to  the  transportation 
system  planning  problem. 

The  subtour  elimination  method  solves  the  classical  assignment  problem 
consisting  of  equations  [1],  (2|,  and  |3J  by  an  efficient  linear  programming  code. 
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Zero-one  solutions  result  (since  the  problem  structure  is  unimodular);  but,  in 
general,  connected  tours  are  not  obtained.  Branch  and  bound  algorithms  have 
been  devised  to  eliminate  subtours,  assuring  satisfaction  of  equation  [4],  Svestka 
and  Huckfeldt  (1973)  report  favorable  computational  experience  with  this  approach 
solving  60-city  problems  with  mean  run  time  of  80  seconds  on  a UNIVAC  1108 
computer.  However,  their  djj  values  were  randomly  generated  and  the  resulting 
matrix  of  values  was  asymmetric  (d^  f d..).  Bellmore  and  Malone  (1971)  have 
shown  that  the  more  practical,  symmetric  case  (d..  = d.j)  is  more  difficult  to 
solve.  In  particular,  the  symmetric  case  requires  an  efficient  technique  to 
eliminate  subtours  of  length  two. 

Little,  et  al.  (1963)  developed  an  efficient  branch-and-bound  algorithm 
for  solving  the  traveling  salesman  problem.  As  discussed  by  Pierce  (1969), 
this  algorithm  may  be  adapted  to  solve  the  multiple  vehicle  problem  with  addi- 
tional constraints  to  reflect  practical  limitations  such  as  vehicle  capacity  and 
time  schedule.  Problems  with  25  cities  may  be  solved  in  approximately  30  sec- 
onds on  an  IBM  7094  computer,  although  run  time  increases  by  a factor  of  10 
for  each  additional  7-10  cities.  This  particular  algorithm  is  useful  only  for  the 


asymmetric  (d_  f d..)  case.  The  transportation  system  problem  assumes 
syipmetric  costs  for  a vehicle  to  travel  between  two  points  (d. ^ = d..).  Approaches 
such  as  that  developed  by  Little,  et  al.  (1963)  are,  in  general,  ineffective  and, 


for  the  most  part,  inappropriate  in  the  case  of  symmetric  distances.  Complex- 


ities which  characterize  the  symmetric  case  arc  stipulated  in  some  detail  by 


Steckhan  (1970). 
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A variety  of  heuristic  procedures  have  also  been  proposed.  While  they 


offer  the  advantage  of  relatively  low  computation  times  for  large  scale  problems, 
they  cannot  guarantee  an  optimal  solution.  Historically,  the  savings  approach 
initiated  by  Clarke  and  Wright  (1964)  has  been  the  most  thoroughly  researched. 

The  fundamental  concept  used  by  the  method  is  that  an  arc  is  added  to  a vehicle 
tour  if  its  inclusion  would  result  in  an  overall  cost  saving.  The  technique  was 
applied  to  single  and  multiple  terminal  delivery  problems  by  Tillman  (1968  and 
1969),  and  embellishments  to  the  basic  approach  were  suggested  by  Tillman 
and  Bering  (1971),  and  Holmes  and  Parker  (1976).  Several  studies — Gaskell 
(1967),  Christofides  and  Eilon  (1969),  and  Webb  (1972) — evaluated  the  perfor- 
mance of  the  savings  approach  and  found  it  to  be  acceptable  in  comparison  to 
other  procedures.  Other,  more  sophisticated  heuristics  have  been  developed — 

Wren  and  Holiday  (1972),  Gillett  and  Miller  (1974),  and  Lin  and  Kernighan 
(1973)  — and  shown  to  produce  good  results  with  reasonable  computation  time, 
even  for  large  problems.  Golden,  Magnonti,  and  Nguyen  (1977)  combined  the 
savings  approach  with  efficient  computational  procedures  and  claim  that  the  re- 
sulting algorithm  is  characterized  by  very  low  run  times.  It  appears,  however, 
that  their  algorithm  terminates  after  quickly  finding  a "local  minimum"  (i.  e. , a 
solution  which  cannot  be  unproved  by  making  a simple,  pairwise  reallocation  of 
points  to  vehicle  routes). 

Two  additional  approaches  have  been  initiated  in  recent  work.  Held  and 

Karp  (1970),  (1971)  have  shown  that  the  traveling  salesman  problem  may  be 

solved  by  formulation  of  the  appropriate  subproblem  which  defines  a minimum 

spanning  tree  which  also  forms  a connected  cycle  including  all  points  in  the 
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network.  Foster  and  ftyan  (1976)  incorporate  characteristics  of  optimal  solu- 
tions to  routing  problems  (as  defined  by  earlier  work)  in  an  efficient  approach 
which  solves  the  vehicle  routing  problem  using  linear  programming  algorithms 
adapted  to  handle  special  cases  which  arise  in  these  problems. 

A more  detailed  evaluation  of  vehicle  routing  problems  is  not  presented 
here  since  they  are  related  to  the  systems  planning  problem  in  a limited  way. 

The  problems  are  similar  in  the  respect  that  both  involve  large  scale  0,  1 integer 
programming  complications.  The  research  discussed  seems  to  be  that  which  is 
most  closely  related  to  the  current  problem.  Excellent  surveys  of  the  work  in 
this  area  have  been  provided  by  Bellmore  and  Nemhauser  (1968);  Pierce  (1969); 
and  Turner,  Ghare,  and  Fourds  (1974).  A rather  comprehensive  bibliography 
of  related  studies  is  included  in  this  report  for  easy  reference. 

Multicommodity  Network  Flows 

Problems  related  to  multicommodity  flows  in  capacitated  networks  also 
encompass  a large  body  of  literature.  Ford  and  Fulkerson  (1958)  initiated  work 
in  this  area  by  developing  an  efficient  algorithm  to  determine  maximal  multi- 
commodity  flows  in  capacitated  networks.  Gomory  and  Hu  (1964)  studied  the 
problem  of  synthesizing  a network  to  satisfy  time  varying  requirements.  Tomlin 
(1966)  formulated  the  problem  of  satisfying  multicommodity  flow  requirements 
in  an  existing,  capacitated  network  at  minimum  total  cost.  He  studied  both  node- 
arc and  arc-chain  formulations  and  indicated  efficient,  large  scale  programming 
techniques  which  might  be  applied  to  solve  problems  of  realistic  size. 
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Cremeans,  Smith,  and  Tyndall  (1970)  reported  an  approach  which  pro- 
vides significant  insight  into  the  proposed  planning  problem.  The  authors  de- 
veloped a minimum  cost  formulation  of  multicommodity  flows  (with  preassigned 
origins  and  destinations)  subject  to  capacity  and  resource  constraints,  and  de- 
livery requirements.  Given  a commercial/dedicated  mode  network  design, 
minimum  cost  flow  plans  could  be  developed  using  this  general  approach.  The 
possibility  of  solving  the  transportation  system  planning  problem  using  this 
general  approach  is  discussed  more  fully  in  Chapter  III.  Cremeans  et  al. 
utilize  an  arc  chain  incidence  matrix  which  defines  the  arcs  included  in  each 
possible  chain  (path)  from  origin  to  destination  for  all  commodity  flows.  The 
arc-chain  approach  offers  the  advantage  that  shipping  time  can  be  calculated  and 
included  in  the  problem  formulation,  while  the  node-arc  formulation  examined  by 
McPherson  and  O'Hara  (1976)  does  not  provide  this  important  capability.  Given 
the  incidence  matrix  of  a particular  network  design,  the  decision  to  be  made  in 
the  optimization  problem  is  simply  how  much  flow  to  prescribe  for  each  chain. 
The  authors  describe  an  efficient  solution  procedure  which  utilizes  a shortest 
path  algorithm  in  a column  generation  scheme.  The  paper  reports  "encouraging" 
results  time-wise  on  problems  involving  150  commodities,  1000  arcs,  and  50 
resources. 

Weigel  and  Cremeans  (1972)  enhanced  the  initial  model  to  allow  more 
realistic  treatment  of  practical  considerations  such  as  vehicles  required  per 
time  period  and  node  capacity  constraints.  A problem  involving  290  arcs,  42 
resources,  and  nine  origin-destination  pairs  was  solved  on  a CDC  6400  computer 
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in  30  minutes.  The  authors  indicate  that  run  time  is  sensitive  to  the  particular 
problem  formulation  used. 


Vehicle  Scheduling 

A variety  of  vehicle  scheduling  problems  have  appeared  in  the  literature. 
Several  types  are  closely  related  to  the  systems  planning  problem  and  are  briefly 
outlined  here  to  record  similarities  and  to  allow  comparison  of  previous  solution 
approaches  with  the  one  selected  for  application  to  the  current  problem. 

Laderman,  Gleiberman,  and  Egan  (1966)  studied  the  basic  problem  of 
allocating  ships  to  (existing)  routes  using  a linear  programming  (LP)  formulation. 
Their  objective  was  to  satisfy  requirements  for  transporting  a single  cargo  con- 
sidering vessel  capacities  and  the  times  at  which  they  would  be  available. 

Conley,  et  al.  (1968)  also  used  an  LP  model  to  minimize  the  total  cost  of  ship- 
ping a single  commodity  through  a rather  large  network.  A more  detailed  prob- 
lem involving  several  commodities,  transshipment  points,  and  two  aircraft  types 
was  studied  by  Gould  (1971).  His  work,  which  was  sponsored  by  the  Military 
Aircraft  Command  (MAC),  developed  an  LP  model  to  minimize  the  costs  of 
moving  aircraft  over  existing  routes  subject  to  operational  constraints.  Benning- 
ton (1970)  was  also  able  to  formulate  an  LP  model  of  a problem  involving  move- 
ment of  military  forces  over  an  existing  network  to  achieve  a desired  schedule 
of  arrivals.  He  employed  a node-chain  formulation  to  describe  flows  and  also 
included  realistic  details  such  as  vehicle  repositioning.  Collectively,  these 
studies  represent  the  capability  of  LP  models  to  prescribe  material  flows  and 
vehicle  schedules.  The  proposed  systems  planning  problem  requires  prescription 
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of  an  optimal  routing  network  in  addition  to  flow  and  time  considerations. 


Several  papers  have  been  addressed  to  this  more  comprehensive  problem. 

Larnpkin  (1966)  solved  a problem  involving  prescription  of  four  aspects 
of  municipal  bus  service:  route  structure  design,  service  frequency,  detailed 
timetables,  and  bus  schedules.  The  author  simplified  the  problem  for  solution, 
dealing  with  different  facets  in  a piecewise  fashion.  He  used  a heuristic  proce- 
dure to  specify  a reasonable  route  structure. 

Bellmore  (1971)  studied  a similar  problem  of  determining  a time  sched- 
ule and  vessel  routing  which  maximizes  a utility  function.  A fleet  of  dissimilar 
tankers  was  employed  to  satisfy  delivery  requirements  for  a single  commodity. 
Delivery  dates  were  constrained  by  predetermined  limits.  The  LP  model  was 
solved  using  Dantzig- Wolfe  decomposition.  A branch  and  bound  algorithm  (of 
the  Land  and  Doig  type)  was  devised  to  obtain  final,  integer  results,  since  a 
fractional  solution  would,  in  general,  be  obtained  from  the  LP  solution. 

A similar  problem  of  selecting  tankers  with  suitable  characteristics 
and  determining  routes  so  that  cargo  is  delivered  within  a specified  time  range 
at  minimum  total  cost  was  researched  by  McKay  (1974).  Employing  a pragmatic 
approach,  he  defined  routes  using  a heuristic  method  (which  was  not  discussed) 
and  solved  the  remaining  portion  as  an  LP  problem  using  a "selective"  rounding 
procedure  to  obtain  integer  results. 

These  studies  have  been  able  to  resolve  network-flow-scheduling  prob- 
lems related  to  the  transportation  system  problem.  However,  none  would  permit 
treatment  of  the  transit  time  of  individual  shipments.  This  is  a major  limitation. 
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since  shipping  time  is  the  prime  measure  of  service  level  afforded  by  a multi- 
model transportation  system.  Additionally,  the  scheduling  problem  posed  is 
more  dynamic  in  nature.  Delivery  times  to  transshipment  points  cannot  be  pre- 
defined within  useful  limits.  Rather,  the  schedules  of  all  vehicles  which  ser- 
vice a transfer  point  should  be  co-ordinated  to  promote  optimal  service  levels. 

One  study  which  addresses  only  the  (dynamic)  vehicle  arrival/departure 
schedule  was  performed  by  Moberly  and  Gorychka  (1970)  at  the  A FIT.  Designed 
to  improve  the  efficiency  of  the  current  Logair  system,  the  study  minimizes 
shipping  time  for  an  established  route  structure.  The  objective  function  consists 
of  terms  of  two  types:  Wjj  (xjj  - y^  ),  and  W[j  (Xjj  - y^  ).  Decision  variables 
x.j  and  y^  represent  vehicle  departure  and  arrival  times,  respectively,  of  flight  i 
at  point  j.  The  factor  applied  to  ground  time,  (x- ■ - y„  f,  weights  the  objec- 
tive function  to  represent  the  amount  of  cargo  shipped  through  point  j on  vehicle  i. 
Analogously,  weighting  factor  Wj.  represents  cargo  transferred  between  flights 
at  point  j.  Summing  over  all  terms  of  this  type,  the  objective  function  results  in 
minimizing  total  weighted  shipping  time.  Their  LI>  model  incorporates  constraints 
which  assure  at  least  the  necessary  minimum  ground  time  and  flight  time  between 
points  as  well  as  relationships  between  flights  involved  in  material  transshipments. 
Numerical  results  show  significant  improvement  over  l97f>  Logair  schedules, 
which  were  designed  by  negotiations  with  the  contractor. 

Vehicle  scheduling  contributes  to  total  system  effectiveness  only  by  allow  - 
ing some  transshipments  to  be  made  expeditiously.  The  ability  to  exploit  'optimum' 
schedules  and  actually  achieve  potential  ship  time  improvements  relies  upon  ma- 
terial handling  capabilities  which  allow  shipments  to  be  transshipod  expeditiously. 
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In  many  cases,  material  handling  requirements  impose  transshipment  delays  on 
the  order  of  several  days.  Therefore  it  appears  that  vehicle  scheduling  repre- 
sents a marginal  opportunity  to  improve  total  system  effectiveness  (particularly 
with  respect  to  the  problems  encountered  by  the  Air  Force  Logistics  Command 
which  were  used  as  a basis  for  identifying  the  transportation  system  planning 
problem).  It  is  therefore  assumed  in  this  study  that  marginal  improvements 
made  possible  by  vehicle  scheduling  may  be  achieved  by  implementing  the 
Moberly/Gorychka  (1976)  once  ’optimal'  vehicle  routes  and  material  flows  have 
been  defined. 

This  assumption  reduces  the  transportation  system  problem  to  two  com- 
ponents. Each  of  these  is  a classical  problem  type  which  has  been  thoroughly 
researched  as  indicated  by  the  literature  review.  However,  the  primary  diffi- 
culty encountered  by  most  existing  approaches  (particularly  those  used  in  vehicle 
routing)  is  that  there  are  no  procedures  available  to  allow  the  necessary  inter- 
action between  the  two  problem  types.  This  interaction  is  necessary  so  that, 
for  example,  vehicle  routes  may  be  redesigned  to  provide  improved  material 
flow  capabilities.  This  study  initiates  research  on  this  topic  and  indicates 
approaches  which  do  provide  such  capability. 
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CHAPTER  III 


A MODEL  FOR  EVALUATING  ALTERNATIVES 
The  model  presented  in  this  chapter  might  be  applied  in  eases  in  which 
management  needs  a means  of  evaluating  a particular  transportation  system 
plan  devised  from  experience  or,  perhaps,  from  a considered  review  of  histor- 
ical trends  in  an  existing  transportation  system.  Certain  operating  systems 
may  be  so  complex  that  it  may  not  be  possible  to  develop  a realistic  model  to 
determine  an  optimal  plan.  In  some  cases,  the  system  complexity  may  impose 
a mathematical  structure  for  which  no  practical  solution  method  is  available. 

In  other  cases,  management  may  like  to  retain  the  decision  making  function 
they  have  performed  historically,  yet  need  a tool  to  provide  an  objective  mea- 
sure of  the  economic  implications  of  their  plan  and  to  assess  the  ship  time  per- 
formance afforded  by  the  overall  transportation  system. 

In  addition  to  providing  measures  of  the  effectiveness  of  a particular 
plan,  the  model  might  be  used  in  the  decision  making  process  to  select  the  pre- 
ferred plan  from  a set  of  feasible  alternatives.  Alternately,  the  model  might 
point  out  certain  weaknesses  of  a proposed  plan  so  that  marginal  improvements 
could  be  identified  and  incorporated  in  the  system  design  stage. 

Finally,  and  perhaps  most  importantly,  the  model  is  actually  a basic  com- 
ponent of  more  complex,  optimizing  models  which  are  discussed  in  later  chapters. 
As  in  the  heuristic  applications,  the  model  is  used  to  evaluate  the'  particular  sys- 
tem designs  in  the  more  advanced  model  forms. 

21 


The  model  is  described  in  the  first  section  of  this  chapter.  Efficient 


solution  procedures  are  discussed  in  detail  in  the  second  section. 

Model  Description 

Given  a network  of  vehicle  tours  traversed  by  components  of  the  de- 
dicated mode(s),  the  transportation  system  planning  problem  reduces  to  a 
multicommodity  material  flow  problem.  In  actuality,  it  is  not  possible  to  state 
a set  of  necessary  and  sufficient  conditions  for  the  dedicated  mode  network 
aside  from  the  system  objective  of  minimizing  total  costs  subject  to  ship  time 
and  operating  constraints.  The  utility  of  a particular  network  design  cannot  be 
evaluated  until  it  is  complete.  There  are  no  criteria,  other  than  the  system 
criteria,  by  which  this  network  can  be  designed.  In  contrast,  in  the  more 
simple,  multiple  vehicle- routing  type  problems,  the  network  of  vehicle  tours 
may  be  defined  with  respect  to  criteria  which  require  one  vehicle  to  visit  each 
point  each  day,  or,  perhaps,  which  provide  sufficient  vehicle  capacity  to  satisfy 
the  needs  of  each  point.  The  transportation  system  problem  has  no  such  criteria 
for  the  vehicle  network.  Rather,  the  primary  design  criteria  is  to  minimize 
total  system  cost.  Since  commercial  modes  are  also  available  for  use,  the  mater- 
ial flow  problem  must  determine  the  manner  in  which  the  dedicated  mode  network 
is  to  be  used,  and,  as  a consequence,  provide  measures  of  system  cost  and  per- 
formance. 

The  material  flow  problem  must  determine  an  optimal  policy  of  using 

both  dedicated  and  commercial  shipping  modes.  In  order  to  satisfy  the  iinpoi 

tant  measure  of  ship  time  performance,  it  is  necessary  for  the  flow  problem  to 
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completely  define  the  paths  from  source,  s,  to  destination,  d,  for  shipments 
made.  These  paths,  which  are  also  specific  to  the  type  of  material,  m,  being 
shipped,  are  called  flow  chains.  Once  the  flow  chain  paths  are  defined,  the 
amount  of  material  to  be  shipped  across  the  flow  chain  must  be  determined. 

Figure  1 presents  the  multicommodity  material  flow  model,  which  is 
very  similar  to  the  one  proposed  by  Cremeans,  et  al.  (1970).  Total  system 
cost  is  composed  of  the  dedicated  mode  cost  plus  that  of  the  commercial  modes 
as  stated  by  the  objective  function,  equation  [8]. 

The  model  is  a linear  program  in  which  the  decision  variables  are  scaled 
to  indicate  the  portion  of  some  msd  shipment  requirement,  f j,  which  is  accom- 
modated by  a specific  flow  chain.  As  stated,  the  model  allows  direct  s to  d ship- 
ments on  either  commercial  air,  or  commercial  truck,  It  is 

assumed  that  a feasible,  albeit  costly,  solution  is  always  possible  through  the 
exclusive  use  of  commercial  air.  Provision  of  the  more  economical,  yet  more 
time  consuming,  commercial  truck  mode  allows  obvious  trade-offs  to  be  consid- 
ered in  system  design.  Additional  commercial  modes  and/or  complicated  ship- 
ping paths  which  use  both  commercial  and  dedicated  modes  could  be  easily  incor- 
porated in  this  type  of  model. 

Equation  (9]  invokes  ship  time  standards  which  are  the  primary  measure 
of  system  performance.  Constraints  of  this  form  could  be  included  for  each  type 
of  material,  for  specific  msd  combinations,  and/or  for  all  shipments  in  an  ox  er- 
all  system  performance  requirement. 

The  total  amount  of  material  transshipped  through  each  point  in  the  net- 
work may  be  constrained  by  material  handling  capability  and  or  available 
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Figure  1. — The  Multieommodity  Material  Flow  Problem 
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warehouse  space.  Equation  [10]  incorporates  this  type  of  practical  restriction 


on  material  flows.  In  fact,  management  may  wish  to  allow  transshipments  at 
only  a few  points  in  the  network  for  administrative  reasons.  It  should  be  noted, 
however,  that  the  opportunity  to  transship  among  vehicles  greatly  enhances  the 
carrying  capacity  of  the  dedicated  mode. 

Commercial  modes  are  assumed  to  have  infinite  capacity,  but  vehicles 
in  the  dedicated  mode  provide  limited  capacity.  The  amount  of  material  shipped 
on  each  arc  traversed  by  a vehicle  in  the  dedicated  mode  must,  therefore,  ob- 
serve this  capacity  restrain,  as  imposed  by  equation  [11].  Vehicles  in  the  net- 
work provided  would  tend  to  be  used  extensively  since  the  flow  problem  incurs 
no  cost  for  use  of  the  dedicated  mode  0 is  assumed). 

Ship  requirements,  f j,  among  all  msd  combinations  are  assumed  to 
be  known  and  must  be  satisfied  by  the  system  plan.  Equation  [12]  invokes  this 
operating  requirement. 

In  seme  applications,  there  may  be  additional  operating  constraints 
which  need  be  imposed.  However,  the  basic  characteristics  of  problems  of  this 
type  are  included  in  this  model.  One  closely  allied  problem  which  may  be  en- 
countered in  practical  situations  is  the  case  in  which  limited  funds  are  provided 
to  operate  a transportation  system  and  the  best  possible  ship  time  performance 
must  be  obtained  within  the  funding  limitation.  This  problem  could  be  solved 
using  the  material  How  model  and  solution  procedure  by  incorporating  two  minor 
modifications.  The  first  would  redefine  the  objective  in  order  to  minimize  I ^ 
(or,  perhaps,  some  weighted  sum  of  the  I"  ).  The  second  would  incorporate 

the  cost  function,  not  as  the  objective,  but  as  an  operating  constraint. 
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As  formulated,  the  flow  problem  is  a large-scale  linear  program.  The 


structure  of,  this  program  is  analyzed  in  the  next  section  and  an  efficient  solution 
procedure  is  described. 


Solution  Procedure 

The  material  flow  problem  is  expressed  in  matrix  notation  in  f igure  2. 
Matrix  dimensions  are  stated  in  Appendix  E.  Two  particular  characteristics  of 
the  model  are  evident  in  this  formulation.  First,  ship  requirements  impose  a 
large  number  of  constraints  in  practical  problems  for  which  there  are  a large 
number  of  points  and/or  material  categories  in  the  system.  Secondly,  there 
are  an  extremely  large  number  of  columns,  each  of  which  represents  one  possible 
flow  chain  or  shipment  route  for  some  msd  combination  in  the  dedicated  mode  net- 
work. Fortunately,  efficient  computational  devices  are  available  to  simplify  each 
of  these  complexities. 

Generalized  Upper  Bounding 

All  of  the  msd  ship  requirement  constraints  may  be  treated  implicitly  by 
the  Generalized  Upper  Bounding  (GUB)  procedure,  which  is  described  in  detail 
by  Lasdon  (1970),  pages  324-340.  The  procedure  incorporates  the  Revised 
Simplex  method  of  solving  a linear  program.  Since  equation  [12]  stipulates  that 
the  sum  of  all  flows  for  each  msd  combination  must  equal  one,  there  must  be  at 
least  one  flow  for  each  combination.  Each  basic  feasible  solution,  therefore, 
has  at  least  one  column  associated  with  each  msd  ship  requirement.  One  basic 
column  for  each  msd  is  defined  as  the  "key  column"  for  the  msd  set  and  other 
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columns,  if  any,  for  that  set  are  expressed  as  a function  of  the  key  column. 

This  approach  allows  the  ship  requirements  to  be  handled  implicitly  so  that  a 
working  basis  consisting  only  of  nonkey,  basic  columns  is  sufficient  for  repre- 
senting a basic  feasible  solution.  Since  the  inverse  of  the  working  basis,  B-^, 
need  be  updated  at  each  iteration  of  the  Revised  Simplex  method,  computation 
time  is  greatly  enhanced  by  using  the  smaller,  working  basis  which  consists  of 
only  (M  +•  I + Q + 1)  rows  and  columns. 

The  GUB  procedure  also  provides  an  efficient  method  to  determine  the 

initial  basic  feasible  solution.  Initially,  Ss  = T , S*  = b.  sa  = b and 

m m i 1 ’ q q’ 

^msd  = 1’0’  A Phase  I procedure  is  needed  to  replace  the  artificial  variables, 
r 

^msd*  This  replacement  is  readily  accomplished,  since  it  was  assumed  that  the 

commercial  air  mode  is  capable  of  delivering  all  shipments  in  such  a way  that 

all  ship  time  standards  are  satisfied.  Since  the  R^sd  are  initially  defined  as 

the  key  variables  (columns)  and  there  are  no  nonkey  members  in  the  working 

basis,  the  replacement  may  be  accomplished  merely  by  redefining  the  key  vari- 
x*A 

ables  to  be  the  Smscj  . The  working  basis  need  not  be  updated  in  this  process, 
conserving  run  time.  Should  commercial  air  not  provide  a feasible  solution  in 
some  application,  the  standard  "Big  M"  method  should  be  applied  in  a routine 
Phase  I process. 

The  basic  feasible  solution  which  results  from  the  suggested  Phase  I 
procedure  is  very  costly,  since  all  shipments  are  accommodated  by  premium 
transportation.  A second  stage  in  the  solution  process  solves  the  material  flow 
problem  using  only  commercial  modes.  The  value  of  this  solution,  BQ, 
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represents  an  upper  bound  on  the  total  system  cost,  since  dedicated  modes  would 
be  specified  only  if  they  can  reduce  total  cost  below 

In  reporting  computational  experience,  Lasdon  (1970)  states  that  the 
Generalized  Upper  Bounding  procedure  may  decrease  computer  run  time  by  a 
factor  of  10  (in  comparison  to  the  straightforward  Revised  Simplex  method).  A 
problem  with  2813  variables,  the  equivalent  of  780  msd  combinations,  and  39 
rows/columns  in  B-^  was  solved  on  an  IBM  7094  computer  in  15  minutes.  Run 
times  of  this  order  of  magnitude  are  certainly  necessary  to  successful  applica- 
tion of  the  material  How  model. 

Column  Generation 

The  large  number  of  columns  associated  with  material  flow  chains-- 
paths  from  s to  d in  the  dedicated  network— may  be  treated  implicitly  by  a column 
generation  procedure.  This  approach  does  not  require  explicit  enumeration  of 
all  possible  flow  chains  in  the  network;  rather,  flow  chains  are  defined  by  solving 
a set  of  shortest-path  subproblems  to  identify  the  column  to  enter  the  basis  at 
each  iteration  of  the  Revised  Simplex  method. 

This  basic  approach  was  first  suggested  by  Ford  and  Fulkerson  (1958). 
Cremeans,  Smith,  and  Tyndall  (1970)  and  Weigel  and  Cremeans  (1972)  report 
successful  application  of  the  approach  in  solving  large-scale,  multicommodity 
network  flow  problems.  The  most  recent  paper  reports  run  time  on  the  order  of 
20  minutes  to  solve  a problem  with  341  rows  on  a CDC  0400  computer.  Since 
the  Generalized  Upper  Bounding  procedure  is  also  being  applied  in  the  current 
research,  the  computational  burden  associated  with  the  large  number  of  rows  in 
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that  example  would  also  be  ameliorated  using  the  approach  discussed  here. 


The  column  generation  procedure  is  described  in  detail  in  the  remainder 
of  this  section.  A typical  column  in  P,  K^,  includes  ship  time,  transshipment 
point  load,  flow  chain  designation,  and  msd  set  membership  information: 

KZ  ~ l (°*  • * Ti!  * fmsd’  * • 0 rz  fmsd)  $1L  fmsd*  * • tIf/  *msd^ 

(alZ  *msd,',aQJt  *msd)  ^1J t •••gmsd’X  ) ^ 

If  the  columns  are  defined  explicitly,  the  t.^  , a ^ , and  gn£  need  be 
specified  for  each  flow  chain: 

ship  time  for  material  type  m on  flow  chain  X 
t^  = 1 if  flow  chain  JL  transships  at  point  i,  0 otherwise 

acj£  = 1 if  flow  chain  X traverses  vehicle  arc  q,  0 otherwise 

g n£  — 1 if  flow  chain  X is  associated  with  msd  = n,  0 otherwise. 

Since  the  columns  are  to  be  generated,  each  of  these  must  be  treated  as  an 
unknown  decision  variable  in  the  subproblems  which  define  flow  chains. 

A t each  iteration  of  the  Revised  Simplex  method,  each  row  has  an  associ- 
ated Simplex  multiplier  calculated  from  B 1 . (CR  is  the  vector  of  objective 
function  coefficients  for  basic  variables  arranged  as  the  columns  of  B. ) The 
vector  of  simplex  multipliers  is  defined  as 


Since  the  objective  is  to  minimize,  the  column,  j,  which  maximi/t-.- 


"(Zj  - Cj)  > 0"  enters  solution  at  each  iteration.  Using  the  formulation  in  l n 
ure  2,  CR  is  a vector  of  zeroes,  except  the  last  element  is  one.  rh<  .-mi;  < 
multipliers  defined  by  CjjB'1  are  therefore  the  elements  in  the  last  r 1 1 

at  each  iteration.  The  column  criterion, 
max  j CB  IT1  - C£|  > 0 
may  be  formulated  using  equations  [14]  and  [15]: 


max 

l 


(Qm  +aM ) * fmsd^  * ^ i ^‘2  *>«sd* 

Q 


msd 


> 0 


+ qfl  ^ (a«l  msd)  + <°>  d0  + £ Bj  «j i ) - 
Ship  time  on  flow  chain  2 may  be  expressed  in  detail  as: 

rP  * °.m  + „1  Tq  »,«.  + Tfm  «lt  * ,4  'll  \l  * «dm 

Since  each  flow  chain  is  related  to  some  specific  msd  combination,  let 
msd 

V * »,«,*> 


1 16] 


1 17  I 


[181 


The  approach  allows  a slack  variable  (Ss,  Sl,  Sa  ) to  enter  into  solution 
if  its  corresponding  Simplex  Multiplier  (om.  Bp  >q  ) is  positive.  All  am,  Bp 
and  yq  are  therefore  negative  in  this  column  generation  procedure.  5^  may  be 
either  positive  or  negative,  although  if  it  is  nonpositive,  column  2 could  not 
enter  solution,  since  it  could  not  "price  out"  at  a positive  value.  Combining 
liquations  [1(1],  |17|,  and  |18],  using  absolute  value  signs,  and  defining 
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I al  I - ! am  I + I aM  I * 
the  criterion  to  determine  the  entering  column  is: 


tiA  « aq £ > X.  ^ = 0, 1 and  form  a connected 

path  from  s to  d.  [19] 

The  value  of  this  criterion  for  each  msd  combination  may  be  found  by 
solving  a shortest  path  problem.  Figure  3 depicts  a typical  subproblem  of  this 
type.  Two  vehicles  are  provided  in  the  dedicated  mode  network.  Vehicle  A 
traverses  the  route  which  visits  points  1,  2,  and  3 in  that  order;  and  vehicle  B 
services  the  route  consisting  of  points  5,  4,  and  2.  Overnight  ship  time  delays 
are  represented  by  the  arcs  between  nodes  1A  and  1A'  as  well  as  between  5B  and 
5B'.  Transshipments  may  occur  at  point  2.  Arc  'distances'  or  costs  are  labeled 
according  to  equation  [19]  components.  Suppose  a shipment  of  material  m begins 
at  source  s = 1 and  is  destined  for  d = 4.  Arcs  from  a 'dummy'  source,  S,  to 
s - 1 and  from  d - 4 to  'dummy'  destination,  D,  account  for  the  constants  re- 
lated to  the  msd  combination.  Dotted  arcs  representing  transshipments  are 
I assigned  'distances'  according  to  the  coefficients  of  tjj,  solid  lines  showing  the 
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arcs  traversed  by  vehicles  are  assigned  the  'distances'  given  by  the  aq^  coeffic- 
ients. The  shortest  path  from  S to  D would  be  the  set  of  arcs  connecting  points 
S-1A'-2A-2B-5B-5B'-4B-D  in  that  order.  The  path  includes  'distance'  con- 
stants associated  with  all  msd  paths,  utilizes  both  vehicles,  incurs  one  over- 
night ship  time  delay  and  one  transshipment  delay.  This  path  is,  in  fact,  the 
only  feasible  one  in  this  example  since  any  other  path  would  include  a subcycle 
(a  visit  to  some  point  more  than  once)  which  would  unnecessarily  increase  ship 
time  and  require  vehicle  capacity  which  might  be  used  for  other  shipments. 

At  each  Revised  Simplex  iteration,  a shortest  path  subproblem  need  be 
solved  for  each  msd  combination  to  determine  the  best  column  to  enter.  Given 
the  path  defined  for  the  msd,  the  column  criterion  is  given  by  equation  [19]. 

The  column  which  gives  the  maximum  positive  criterion  i.  the  one  to  enter  at 
this  iteration.  If  a column  has  a negative  criterion  value,  it  is  not  a candidate 
to  enter;  and  if  the  shortest  paths  for  all  msd  yield  nonpositive  criteria,  the  cur- 
rent solution  is  optimal  since  no  column  could  enter  and  improve  the  current 
solution.  A column  already  in  solution  would  not  be  generated  again  by  this  pro- 
cedure since  (by  definition  of  the  Simplex  Multipliers)  it  would  'price  out'  at 
zero  and  would  therefore  not  be  a candidate  to  enter. 

In  lieu  of  solving  msd  shortest  path  problems,  a single  subproblem  ap- 
propriately defined  could  be  used  to  solve  2,  3, ... , msd  of  the  original  subprob- 
lems simultaneously.  Figure  4 depicts  a network  in  which  all  msd  subproblems 
could  be  solved  as  a single  shortest  path  problem.  A mathematical  statement  of 
this  single  problem  is: 
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AGGREGATED  SHORTEST  PATH  NETWORK 
FOR  COLUMN  GENERATION 


Min  Z 


Subject  to: 

• I 

2 S bl  ~ Cl  ~ fmsd[  ( I I Tim  + I Pi  I > hi 

Q _ t 

+ S ( I I Tq  + I >q  I \l 

q~i 

i 

+ .=  <l°l  I 4.)  X,A 
1^1 

+ I I < °sm  + Udm  ) 

Z - 0,  tj£,  , X.^  - 0, 1 and  form  a connected  path 

from  s to  d. 

Computational  experience  is  necessary  to  determine  the  best  number  of  subprob- 
lems to  solve  simultaneously  in  this  formulation. 

In  any  event,  the  usefulness  of  this  overall  solution  approach  is  largely 
dependent  upon  solving  the  subproblem(s)  efficiently.  A special  algorithm  was 
devised  to  conserve  computer  storage  space,  and  to  provide  efficient  computation 
times.  The  algorithm  is  based  on  an  early  paper  by  Dijkstra  (1959)  in  which  he 
suggested  an  approach  which  has  apparently  (Dreyfus  (1967),  Golden  (1976), 
Elmaghraby  (1970)  ) not  been  bettered  by  more  recent  research  on  the  shortest 
path  problem. 

The  algorithm  solves  a set,  p,  of  the  msd  subproblems  simultaneously 
by  discovering  the  shortest  path  from  S to  D as  shown  in  Figure  4.  Processing 
steps  operate  on  the  nodes  N in  the  set  of  subnetworks,  placing  them  on  list  A if 
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the  shortest  path  from  S to  N has  been  discovered,  on  list  B if  some  path  from 


S to  N has  been  defined,  and  on  list  C if  no  path  from  S to  N has  yet  been  defined. 

To  begin,  all  nodes  for  msd  C r|  are  on  list  C.  The  source  nodes,  s,  are  removed 
and  assigned  to  list  A and  all  nodes  accessable  from  the  source  nodes  are  placed 
on  list  B.  Entries  on  list  B are  ordered  low  to  high  with  respect  to  the  distance 
from  S to  N,  At  each  iteration,  the  first  entry  on  list  B,  N,  is  removed  and 
placed  on  list  A since  the  shortest  path  from  S to  N is  now  evident — any  other 
path  from  S to  N would  necessarily  be  longer.  Each  node  J which  is  accessable 
from  node  N is  then  checked  to  see  if  the  path  through  N to  J is  better  than  any 
identified  previously.  List  B and,  if  necessary,  list  C are  updated  accordingly. 

The  shortest  path  from  S to  D is  known  when  any  of  the  d ^ for  msd  C r) 
reaches  the  top  of  list  B.  The  column  criterion  may  then  be  calculated  to  deter- 
mine if  this  flow  chain  should  become  the  next  entering  column.  If  no  flow  chain 
which  satisfied  equation  [19]  is  defined  in  this  procedure,  the  current  solution  is 
optimal  for  the  material  flow  problem. 

In  addition  to  using  the  Dijkstra  approach  as  the  underlying  mechanism 
in  the  shortest  path  algorithm,  several  additional  features  are  incorporated  to 
promote  efficiency.  Arc  'distances’  need  not  be  calculated  until  needed,  so  the 
'distances'  of  some  number  of  arcs  may  never  be  calculated.  Nodes  and  arcs 
need  not  be  added  to  represent  transshipments  and  overnight  vehicle  terminal 
delays  as  suggested  by  the  conceptual  presentation  in  Figure  3.  Rather,  these 

factors  may  be  treated  without  expanding  network  size.  Finally,  an  interval  bi- 

•* 

section  procedure  may  lie  used  to  reorder  nodes  on  list  B.  Tills  procedure  is 
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expected  to  be  efficient  in  solving  large  problems.  A detailed  statement  of  the 
algorithm  follows. 


Step  0:  Initialization 

A.  Initially,  all  nodes,  I,  for  msd  combinations  C q are  on  list  C with 
D(I)  = 0.  A vector  of  data  for  each  vehicle  arc  defines  the  beginning 
node,  I;  the  ending  node,  J;  the  associated  vehicle,  k;  the  number  of 
the  successor  arc  traversed  by  vehicle  k,  NSUC;  the  number  of  trans- 
shipments possible  at  point  J and  the  successor  vehicle  arc  numbers 
for  each  transshipment. 


Step  1:  Originating  Flow  Chains 


A.  For  each  msd  C q 

Remove  the  source  node,  s,  from  list  C.  Calculate 


D(s)  ~ ~ CJ l + fmsd  I amsd  I °sm 

If  D(s)  s 0,  drop  this  msd  combination  since  it  could  not  provide 

the  entering  column.  Otherwise,  add  node  s to  list  A,  since  D(s)  is 
the  shortest  length  path  from  S to  s.  Process  each  arc  q which 
emanates  from  node  s in  subnetwork  msd  according  to  step  1/B. 

B.  For  each  arc  emanating  from  s in  msd  C q 

Calculate  the  distance  along  arc  q and  at  its  end  point,  J, 

Bq  ' < I "rnsd  I Tq  + I % I I 

If  point  J is  the  destination,  dmscj.  add  | omg(j  | U to  B^  to  account 
for  unload  time  at  the  destination.  If  J / (l,nS(p  and  J is  the  termination 
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point  for  the  k tour,  add  | c? mSci  | to  to  account  for  overnight 

ship  time  delay  before  next-day  movement.  Calculate  distance  measure 
D(J)  ^ D(S)  + fmsd  Dq. 

If  D(J)  — 0,  drop  this  potential  path,  since  it  would  not  be  able  to  pro- 
vide the  shortest  (nonpositive)  path  from  S to  J.  Otherwise,  remove  J 
from  list  C and  place  it  on  list  B in  order  of  increasing  D(.  ) for  the 
nodes  on  list  B,  Record  the  msd  subnetwork  associated  with  this  node 
(NET)  and  the  predecessor  of  node  J on  list  A,  PRED(J)  - s. 

C.  After  all  source  nodes  for  msd  C q and  the  vehicle  arcs  which  emanate 
from  them  have  been  processed,  go  to  step  2. 


Step  2;  General  Iterative  Process 

A.  If  list  B is  empty,  stop;  the  entering  column  cannot  be  from  the  set  q. 
Otherwise,  remove  the  first  entry  from  list  B,  node  N.  The  shortest 
path  from  S to  N is  now  known  to  be  the  stored  value,  D(N),  the  prede- 
cessor node  on  this  path  is  PRED(N),  the  applicable  msd  subnetwork  is 
the  stored  value,  NET.  Add  node  N to  list  A.  If  N is  the  destination, 
dm  j,  go  to  step  4.  If  there  is  no  successor  arc,  NSUC,  for  the  vehicle 
which  brought  the  flow  chain  to  N,  go  to  step  3.  Otherwise,  go  to 
step  2/B. 


B.  The  successor  arc  q for  vehicle  k begins  at  point  N and  ends  at  point  J. 
Determine  the  ’distance'  along  arc  q and  at  its  end  point: 

i>„  < I "msd  I T‘  ' I y„  I > • 


As  in  step  l/B,  update  Dq  if  point  J is  the  destination,  dmgd,  or,  if 
J f dmgd  and  J is  the  termination  point  for  vehicle  k. 

Calculate  distance  measure 
D(J)  = D(N)  + Dq. 

If  D(J)  - 0,  drop  this  potential  path  to  J and  go  to  step  3. 

C.  If  node  J in  this  msd  subnetwork  is  on  list  C,  remove  it  and  order  it  on 
list  B using  measure  D(J).  Record  the  msd  network  (NET)  for  J and 
set  PRED(J)  = N.  Go  to  step  3. 

D.  If  node  J in  this  msd  subnetwork  is  already  on  list  B and  if  D(J)  cal- 
culated above  is  greater  than  the  stored  distance  from  S to  J,  D(J),  go 
to  step  3.  Otherwise,  a shorter  path  from  S to  J has  been  found. 

Give  D(J)  its  new  value  calculated  above,  reorder  J on  list  B,  and 
record  PRED(J)  = N.  Go  to  step  3. 

3:  T rans  shipments 

A.  If  there  can  be  no  transshipments  at  node  J,  return  to  step  2 /A. 
Otherwise  process  each  of  the  vehicle  arcs  which  transship  out  of 
node  J according  to  the  following  substeps. 

B.  The  flow  chain  leaves  vehicle  k at  point  J,  incurs  delay  Tjm>  and 
departs  J on  vehicle  k'  on  arc  q bound  for  point  J'. 

Determine  the  distance  along  arc  q and  at  its  end  point: 

“q  1 I "t  1 1 Jm  + Pj  1 + < I "l  I Tq  ' K I ' 

Update  D(|  lor  II and  T ^ (if  appropriate)  as  in  step  2/B. 


r 

I II 

Calculate  distance  measure 
D(J’)  = D(J)  + Dq. 

If  D(J')  - 0,  drop  this  potential  path  to  J'  . 

Otherwise,  process  node  J'  using  D(J')  according  to  steps  2/C  and  2/D. 

C.  Return  to  step  2/A  after  all  possible  transshipments  out  of  J have  been 
examined. 

Step  4;  Destination  Node  dmgd  has  been  reached. 

A.  Node  dmgcI  was  just  removed  from  list  B in  step  2/A,  the  shortest 
path  from  S to  any  destination  for  the  msd  __  r\  has  thus  been  identi- 
fied (the  algorithm  generates  paths  from  S to  nodes  in  increasing 
distance  so  no  path  to  another  destination  could  be  better  than  the  one 
just  discovered).  Trace  the  path  from  d d back  to  S using  the 
PRED(N)  data.  Accumulate  the  total  ship  time  along  this  flow  chain 
as  defined  by  equation  [17),  and  also  accumulate  the  total  column 
criterion  as  defined  by  equation  [19].  If  the  column  criterion  is  greater 
than  the  current  best,  record  the  new  column  data  as  defined  in 
equation  [19). 

B.  If  all  msd  have  been  included  in  some  subproblem,  stop.  (Either  the 
best  entering  column  has  been  defined  or  it  has  been  shown  that  no 
column  will  enter  and  the  current  solution  is  optimal).  Otherwise, 
increment  the  set  q to  include  a new  set  of  the  msd  combinations  and 
go  to  step  0. 
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CHAPTER  IV 


AN  OPTIMIZING  APPROACH 

A hypothesized  economic  relationship  between  the  cost  of  the  dedicated 
mode  vehicle  network  and  the  total  system  cost  is  depicted  in  Figure  5.  If  no 
dedicated  mode  be  provided,  the  minimum  cost  shipping  plan  using  commercial 
modes  exclusively  is  BQ  as  described  in  the  previous  chapter.  Provision  of  a 
'small'  dedicated  mode  network  may  actually  increase  total  system  expenses, 
since  vehicles  in  the  dedicated  mode  may  not  be  used  to  capacity.  As  the  scale 
of  the  dedicated  mode  is  increased,  more  transshipments  become  possible  and 
material  flow  is  facilitated.  A variety  of  network  designs  may  be  possible  for 
any  particular  level  of  expense.  However,  each  design  may  give  rise  to  differ- 
ent material  flow  characteristics  and,  hence,  to  different  system  cost.  This 
phenomenon  is  reflected  by  the  band  in  which  total  costs  are  hypothesized  to 
lie.  As  the  dedicated  mode  is  made  larger  and  larger,  capacity  is  available  to 
transport  most  shipments  and  transshipment  opportunities  are  present  to  pro- 
mote smooth  flow.  Continued  increases  in  scale  at  that  point  bring  about  little, 
if  any,  system  improvement;  and,  ultimately,  excess  capacity  beyond  that 
needed  is  indicated.  It  is  expected  that  a rather  broad  range  of  network  designs 
would  provide  total  system  cost  near  the  optimum.  The  relative  locations  of 
points  to  be  served  in  combination  with  other  important  features  evident  in  a 
particular  application  may  allow  an  analyst  to  determine  a solution  within  this 
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broad  range  using  only  the  heuristic  procedure  described  in  the  previous 
chapter. 

However,  in  other  applications  it  may  be  difficult  to  design  a good  net- 
work with  heuristics  and  a method  to  determine  an  optimal  solution  is  necessary. 
Synthesizing  a network  design  is  a complex  decision  process.  Vehicle  tours 
must  be  determined  by  material  flow  requirements  and  total  system  cost.  Char- 
acteristics of  the  less  complex,  multiple- vehicle  routing  problems — providing 
one  vehicle  to  service  each  point  daily,  designing  the  network  so  that  the  total 
travel  distance  is  minimized,  using  cloverleaf  route  geometries — may  have 
little  bearing  on  practical  problems  of  the  transportation  system  planning  type. 

A mathematical  model  which  does  specify  an  optimal  transportation 
system  plan  is  presented  in  this  chapter.  All  design  features  described  in  the 
problem  statement  are  incorporated  in  the  model,  which  is  a large  scale,  mixed 
0, 1 integer,  linear  program. 

The  model  is  presented  in  the  first  section  of  this  chapter.  A solution 
algorithm  designed  to  exploit  the  particular  structure  of  the  model  is  described 
in  the  second  section. 

Model  Description 

The  model  requires  feasible  vehicle  tours  to  be  defined  and  made  avail- 
able for  use  in  designing  the  network.  The  0,  1 integer  portion  of  the  problem 
then  determines  which  set  of  vehicle  tours  to  use  in  the  optimal  network  design. 
The  optimality  of  the  resulting  system  plan  is,  of  course,  relative  to  the  set  of 
tours  provided.  However,  this  may  not  be  particularly  restrictive,  since 
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experience  should  be  a good  indicator  of  the  types  of  vehicle  tours  most  likely  to 


service  system  needs  efficiently  and  since  a large  number  of  tours  could  be 
made  available  for  possible  use  by  the  design  process.  For  problems  with  a 
small  number  of  points,  all  feasible  tours  could  be  enumerated  using  a simple 
algorithm  described  later.  At  worst,  this  approach  allows  the  analyst  to  eval- 
uate a set  of  design  combinations  far  in  excess  of  the  number  possible  through 
heuristic  or  manual  procedures. 

A detailed  statement  of  the  model  appears  in  Figure  6.  The  objective 
function  combines  the  cost  of  material  flow  and  the  expense  of  vehicle  tours  used 
in  the  dedicated  mode.  Tour  cost,  C may  reflect  a number  of  components 
including,  for  example,  a fixed  investment-type  expense,  a cost  to  travel  be- 
tween each  of  the  points  on  the  tour,  and  costs  to  service  each  point.  Equations 
[21],  [22],  and  [24],  which  provide  constraints  on  ship  time,  transshipment  point 
loading,  and  ship  requirements  are  the  same  as  equations  [9],  [10],  and  [12]  in 
the  material  flow  model  of  Figure  1. 

The  remainder  of  the  model  accounts  for  the  network  design  process. 
Equation  [23]  defines  vehicle  arc  capacity,  but  differs  from  equation  [11]  in  the 
material  flow  model  in  one  important  respect.  It  is  assumed  in  this  model  that 
each  vehicle  provided  establishes  a set  of  constraints  of  the  form  in  equation  [23]. 
Arcs  not  used  in  a particular  network  design  cannot  be  used  for  material  flow 
and  have  right  hand  side  equal  to  zero.  Each  arc,  q,  is  defined  by  a vector  of 
information  including  r<  beginning  and  ending  points,  the  specific  vehicle  in- 
volved, and  whetln  • • tlx  'n  <t  u second  time  the  vehicle  travels  from  the 
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Complementary  Vehicle  Relationships 
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Vehicle  Network  Constraints 

As  required  in  a specific  application  to  restrict,  for  example, 
the  total  vehicle  capacity  provided,  the  total  number  of  vehicle 
visits  to  a point,  or  to  require  a return  tour  supplementing 
vehicle  paths. 


Non-negativity  and  Integer  Requirements 
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L ' msd’  msd 


= 0,1 


Vi, , msd 


V k,  rCT(k) 


Figure  6.  — Continued 


beginning  to  ending  point.  This  level  of  detail  is  necessary  to  assure  that  no 


vehicle  is  overloaded  on  any  arc. 

The  remaining  constraints  relate  to  vehicle  tours.  Equation  [25]  re- 
quires each  vehicle  to  be  used  in  at  most  one  tour  in  a network  design.  Should 
the  null  tour  (the  vehicle  actually  does  nothing)  be  defined  for  a vehicle,  equa- 
tion [25]  should  be  stated  as  a strict  equality,  indicating  that  each  vehicle  must 
have  some  tour  in  solution. 

Equation  (26[  provides  a means  of  improving  solution  times  in  some 
cases  by  defining  vehicles  in  predetermined  pairs  to  be  complementary  to  each 
other.  This  device  permits  several  possibilities.  First,  the  null  tour  for 
vehicle  k could  be  formulated  as  the  only  tour  allowed  for  k,  providing  a means 
of  invoking  the  relationship  noted  above.  Secondly,  this  type  of  constraint  could 
distinguish  between  tour  paths  and  tour  routes  traversed  by  a single  vehicle. 

For  example,  vehicle  k may  originate  at  Los  Angeles,  visit  several  points  and 
terminate  at  New  York.  The  set  of  all  such  paths  could  be  viewed  as  the  possible 
tours  for  k.  An  alternate  network  strategy  might  call  for  k to  traverse  a route, 
traveling  eastward,  say,  to  Dallas  and  subsequently  terminating  at  its  point  of 
origin,  Los  Angeles,  to  complete  the  route.  In  actuality,  there  may  be  only  one 
vehicle  and  the  complementarily  constraint  indicates  that  only  one  type  of  tour  is 
feasible.  Additionally,  the  constraint  may  improve  solution  time  by  providing  a 
logical  means  of  branching. 

Should  a tour  path  be  defined,  other  tour  paths  also  need  be  defined  to 
return  the  vehicle  from  its  termination  point  the  next  day.  Otherwise,  one  point 
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may  be  a 'sink'  at  which  all  vehicles  terminate  and  therefore  originate  the  next 
day.  This  and  other  constraints  which  management  might  impose  on  network 
design  in  any  particular  application  are  represented  by  the  note,  equation  |27], 

Each  of  the  defined  tours  is  assumed  to  satisfy  all  constraints  imposed 
by  management  as  well  as  any  limitations  which  result  from  vehicle  operating 
characteristics.  Examples  of  such  constraints  are:  the  number  of  points  visited 
by  a vehicle;  the  tour  length  in  time  and/or  distance  stipulated  by  vehicle  capa- 
bilities, or,  perhaps,  crew  requirements;  and  the  number  of  times  a vehicle 
visits  a particular  point.  Furthermore,  it  is  assumed  that  each  vehicle  tra- 
verses its  assigned  tour  during  the  course  of  a day.  Routes  are  traversed  daily 
and  in  the  case  of  paths,  vehicles  are  redefined  and  depart  from  their  termina- 
tion point  the  next  day.  Alternate  assumptions  concerning  tour  duration  and 
frequency  could  be  incorporated  in  the  model  if  necessary. 

Even  though  the  model  may  be  stated  using  only  a few  types  of  constraints, 
it  does  represent  a large  scale  mathematical  programming  problem.  Available 
mixed  0, 1 integer,  linear  programming  codes  might  be  applied  to  solve  problems 
on  limited  scope.  But  this  approach  is  expected  to  find  limited  success  unless 
the  GUB  and  column  generation  procedures  described  in  the  previous  chapter 
can  be  applied.  An  alternate  procedure  which  does  allow  use  of  these  efficient 
techniques  results  from  application  of  the  concepts  of  Benders'  decomposition. 
These  concepts  were  originally  developed  by  Benders  (1962)  and  were  later  ex- 
tended by  Geoffrion  (1972).  The  technique  is  well  summarized  by  Lasdun  (19701. 
Benders'  decomposition  allows  the  system  model  in  Figure  6 to  be  re-expressed 
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as  two  interacting  models.  One  of  these  problems  is  defined  as  Problem  LP 
in  Figure  7;  and  the  other,  Problem  01,  is  stated  in  Figure  8.  Each  of  these 
problems  is  more  simple  than  the  original  problem,  and  efficient  solution  pro- 
cedures may  be  devised  to  exploit  the  special  structure  of  each  providing  effic- 
ient computation  capability. 

Problem  LP  is  essentially  the  same  as  the  multicommodity  material 
flow  problem  in  Figure  1.  LP  accepts  a vehicle  network  (specified  by  01)  and 
determines  the  least-cost  shipping  plan  using  both  dedicated  and  commercial 
modes.  In  turn,  LP  provides  information  to  01  by  adding  a constraint  which 
essentially  gives  an  indication  of  the  relative  utility  of  the  vehicle  network.  The 
solution  procedure  continues  iteratively  in  this  fashion  until  the  optimal  solution 
is  determined.  Convergence  is  guaranteed  by  the  theory  of  Benders'  decompos- 
ition, but  a large  number  of  iterations  may  be  required.  However,  the  litera- 
ture records  successful  application  of  the  method  to  a variety  of  problems  (see, 
for  example,  Geoffrion  and  Graves  (1974)  and  Richardson  (1976)  ).  Considering 
these  reports  of  computational  experience,  the  decomposition  approach  is  ex- 
pected to  provide  an  efficient  solution  procedure  for  the  transportation  system 
planning  problem. 

The  right-hand  sides  of  constraints  represented  by  equation  [32]  are  zero 
if  the  arc  is  not  used  in  a given  vehicle  network  and  equal  to  bq  if  arc  q is,  in 
fact,  defined  in  the  network.  At  each  /c-iteration,  equation  [32]  need  be  defined 
only  for  those  arcs  in  the  vehicle  network  defined  by  Problem  01.  Zero  arc 
capacity  merely  indicates  that  no  material  How  chain  may  use  the  arc.  Excluding 
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Min 


p = f 

o msd 


' L 
£ Ct 

i-1  * 


pf>  + 


v crA  SrA  + 
msd  msd  msd  msd 


v rrT  rT 
msd  msd 


1291 


(30) 


[31) 


132] 


[331 


i>(t) 

l 


, rA 

''’msd 


rT 

‘ msd 


’ 0 


Vi. , msd 


1 34 1 


l i^ure  7.  — Problem  1,1’ 
31 


it  from  the  problem  obviously  accomplishes  the  same  result.  As  shown  in 


Appendix  F,  the  dual  variables  associated  with  these  undefined  arcs  must  equal 
zero  whether  or  not  they  are  included  in  the  formulation  of  LP. 

The  objective  function  of  the  dual  of  LP  is 

M ! _ Q 

max  2 r a + 2 b.p.  + 2 b 

m=l  m m i-1  1 1 q=l  q 

In  most  applications,  the  dual  of  the  linear  programming  problem  would  be 
solved  rather  than  the  primal,  since  the  set  of  feasible  solutions  in  the  dual  is 
not  affected  by  the  values  assigned  to  the  binary  variables.  The  advantage  of 
using  the  dual  is  that  the  solution  found  at  the  last  k -iteration  may  be  used  as 
the  initial  basic  feasible  solution  for  the  current /c-iteration.  Only  the  dual  ob- 
jective function  coefficients  change  as  a function  of  the  binary  variables  as  noted 
by  equation  [35].  However,  in  the  current  application,  it  is  recommended  that 
the  primal,  LP,  be  used  to  take  advantage  of  several  possibilities  to  reduce 
computation  time:  1)  only  constraints  with  nonzero  right-hand  sides  need  be 
defined  in  equation  [32]  at  each  iteration,  2)  the  GUB  method  may  be  applied  to 
the  primal,  3)  the  column  generation  scheme  described  in  the  previous  chapter 
may  be  applied  to  the  primal  problem.  Since  an  efficient  Phase  I procedure  is 
available,  use  of  the  primal  appears  preferable  to  use  of  the  dual  in  this  application. 

It  is  well  known  that  the  optimal  solution  to  LP  has  the  same  value  as 
the  optimal  solution  to  its  dual  problem.  Thus,  for  optimal  values  of  all 
variables: 
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Considering  the  directions  of  the  inequalities  in  the  constraints  of  problem  PI, 
it  is  seen  that  a , P , and  y are  nonpositive  and  & is  unrestricted  in  sign.  The 
value  of  the  dual  objective  function  is,  however,  nonnegative  as  noted  by  equa- 
tion [36]. 

At  each  /c-iteration,  the  dual  objective  function  is  combined  with  the  ex- 
pression for  vehicle  tour  cost  to  form  an  additional  Benders'  cut  (or  constraint) 
in  Problem  01: 
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equation  [371  may  bo  re-expressed  as  in  equation  |43|. 
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Problem  01 


The  other  constraints  in  Problem  01  are  taken  from  Problem  PI  and 


carry  the  same  implications  as  discussed  previously.  The  continuous  variable, 

Z,  is  defined  by  the  Renders'  cut  which  places  the  highest  cost  assessment  on  a 
particular  network  design.  Benders'  cuts  essentially  place  a value  on  a proposed 
network  which  is  determined  from  network  designs  previously  evaluated  by 
Problem  LP.  In  particular,  y = 0 if  the  vehicle  which  traverses  arc  q is  not 
loaded  to  capacity.  The  cuts,  therefore,  tend  to  feed  back  more  information  to 
01  if  a number  of  vehicle  arc 5 are  loaded  to  capacity.  Fortunately,  Problem  LP 
attempts  to  load  arcs  to  capacity  since  they  are  available  at  zero  cost,  while  ex- 
penses are  incurred  for  use  of  commercial  modes.  The  cuts  do  not  provide  an 
indication  of  individual  tours  which  might  be  effective.  Rather,  they  place  a 
value  on  the  network,  or  complete  set  of  vehicle  tours.  The  effective  cost  of 
including  an  arc,  which  was  loaded  to  capacity  in  a previous  solution,  is  reduced 
by  the  y terms.  The  y indicate  the  utility  of  a vehicle  arc  with  respect  to  the 

CJ/C  CJK. 

material  flow  subproblem  and  tend  to  identify  arcs  which  may  be  used  effectively 
by  material  flow  if  they  are  provided  by  solution  of  the  network  subproblem. 

In  a general  formulation  involving  Benders'  decomposition,  a second  type 
of  constraint  may  have  to  be  augmented  to  the  integer  problem.  The  function  of 
the  second  type  of  cut  is  to  assure  that  solutions  provided  by  the  integer  problem 
allow  the  primal  linear  program  to  attain  a feasible  solution.  These  cuts  are  not 
required  in  the  current  formulation,  since  it  is  assumed  that  exclusive  use  of 
commercial  air  is  always  feasible.  Thus  LP  may  attain  a feasible  solution  no 
matter  what  vehicle  network  is  specified. 


Solution  Procedure 


Efficient  procedures  to  solve  Problem  LP  were  described  in  Chapter  IEL 
Steps  required  to  solve  Problem  01  and  the  combination  of  01  and  LP  are  de- 
tailed in  this  section. 

Tour  Provision 

A set  of  feasible  tours  for  each  vehicle  may  be  defined  by  the  analyst 
and  made  available  for  use  in  Problem  01.  This  process  would  not  be  expected 
to  be  particularly  difficult  or  time  consuming,  since  tours  could  be  readily  de- 
fined by  picking  likely  candidates,  considering  the  relative  locations  of  the 
points  to  be  served  and  the  material  flow  requirements  between  them.  Tours 
need  not  be  confined  to  limitations  typically  imposed  in  the  multiple- vehicle 
routing  problem.  For  example,  a given  vehicle  tour  might  visit  one  point  a 
number  of  times,  or  all  vehicle  tours  may  service  one  of  the  points  to  assure 
transshipment  capability. 

All  feasible  tours  may  be  enumerated  using  the  following  algorithm: 

Step  0:  Initialization 

Define  the  set  of  0-length  tours  to  be  empty.  Set  n - 0.  Go  to  step  1. 

Step  1:  Tour  Definition 

Define  (n  + 1)- length  tours  by  augmenting  the  termination  point,  t, 

(for  vehicle  k)  to  each  of  the  n-length  partial  tours  which  do  not  have  t 
as  the  last  point.  Check  tour  feasibility  with  respect  to  time  and/or 
distance  length,  number  of  visits/point,  etc.  Record  any  feasible 
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(n  + l)-length  tours.  If  (n  + 1)  = maximum  feasible  tour  length,  stop. 
Otherwise,  go  to  step  2. 


Step  2:  Partial  Tour  Extension 

Increment  n.  Define  n-length  partial  tours  by  augmenting  point  i 
(1=1,  2, ... , I;  i/last  point  on  (n-l)-length  tour)  to  each  of  the  (n-l)-length 
partial  tours.  Check  tour  feasibility  with  respect  to  time  and/or  distance 
length,  number  of  visits/point,  etc.  Drop  any  infeasible,  partial  tours 
from  the  list  of  n-length  partial  tours.  Go  to  step  1. 

Since  this  algorithm  would  define  a large  number  of  tours,  it  should  be  used  only 
when  each  vehicle  can  service  a low  number  of  points  and/or  when  feasibility  con- 
straints rule  out  a large  number  of  potential  tours.  Of  course,  the  algorithm  could 
be  used  to  define  a set  of  tours  from  which  the  analyst  may  pick  likely  candidates 
to  make  available  to  Problem  01. 

It  is  also  recommended  that  the  analyst  define  a 'good'  vehicle  network 
so  that  a relatively  tight  upper  bound  on  the  solution  to  PI  may  be  defined  early. 
This  should  not  be  difficult  assuming  that  the  analyst  has  some  familiarity  with 
historical  operations  of  the  system  under  study  or  of  similar  systems. 

Solution  Bounds 

Iterations  may  begin  by  solving  LP  using  the  initial  network  defined  by 
the  analyst.  Useful  bounds  on  the  solution  to  PI  and  on  the  value  of  / may  be 
defined  at  each  iteration  according  to  the  theory  of  Benders'  decomposition 
(l.asdon,  l!)7(>). 

57 


1 

/c  /c 

If  at  iteration  fc,  (Z  , T ) solves  Problem  01  and  Pq  is  the  optimal 
value  of  Problem  LP  given  T*  , then  [T K , PK  (SrA)'c  , (SrT)'c  ] is  a feasible 
solution  to  the  original  Problem  PI.  The  corresponding  objective  function  value 
in  PI  is 

F*  = Z C , T*  + P K . 
gr  rk  rk  o 

This  relationship  maybe  used  to  define  the  upper  bound  on  the  value  of  Z for  the 
next  iteration: 

UB^1  = min  j UB^ , FK  } . [45] 

A lower  bound  on  the  value  of  Z for  the  next  iteration  is 

LB K+1  = ZK 

since  the  value  of  Z is  monotonically  nondecreasing  for  successive  < -iterations. 

Therefore,  at  the  next  iteration 

LB  'C+1  s zK+1  s ubk+1, 

and  at  the  optimal  solution 

LB K = ZK  = UB K . 

The  iterative  procedure  may  be  stopped  when  [ TK , p\  (SrA)K  , (S1**)*"  ] 

provide  a solution  which  is  within  e of  optimality;  i.  e. , stop  when 

(UB*  * 1 - LB  * 1 1 ) s f . The  true  optimal  solution  is  then  known  to  be 
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somewhere  between  the  lower  and  upper  bounds.  In  problems  which  converge 
slowly,  it  may  be  prohibitively  expensive  to  determine  the  optimum  exactly. 
The  ability  to  define  e optimality  may  therefore  be  a useful  characteristic  of 
this  approach. 


An  Algorithm  for  Problem  01 

Since  n tours  are  provided  for  vehicle,  k,  there  are  a total  of 
K 

K 

N = 2 n. 

k=l  k 

binary  variables  in  addition  to  the  one  continuous  varible,  Z,  in  Problem  01. 

The  total  number  of  combinations  of  N binary  variables  is  2‘s.  However,  since 
each  vehicle  may  use  only  one  tour  at  level  1 in  a feasible  solution,  there  are 
only 

_ K 
N = FI  n, 

k=l  k 

feasible  combinations.  If  complementary  vehicle  relationships  and  other  net- 
work constraints  are  imposed,  the  number  of  feasible  solutions  is  less  than  N. 
Standard  integer  programming  codes  could  be  applied  to  solve  Problem  01  but 
they  would  be  forced  to  operate  with  the  N solution  space  rather  than  that  of  N. 
Additional  binary  variables  would  be  required  to  re-express  Z as  a function  of 
integer  variables,  further  complicating  the  problem.  A special  purpose,  im- 
plicit enumeration  algorithm  was  developed  to  exploit  the  structure  of  Problem  01 
and  provide  an  efficient  solution  procedure.  Since  Problem  01  need  be  solved 
several  times  in  the  process  of  solving  Problem  PI,  an  efficient  solution 


r 

procedure  contributes  to  the  overall  capability  of  the  decomposition  approach  to 
solve  transportation  system  planning  problems. 

The  special  purpose  algorithm  operates  in  a manner  similar  to  that  of 
Balas's  implicit  enumeration  strategy.  Balas's  (1965)  strategy,  in  combination 
with  the  enumeration  logic  devised  by  Geoffrion  (1966)  (see  Taha  (1971)  for  a 
detailed  description)  forms  the  basis  -or  a number  of  efficient  integer  program- 
ming algorithms  (see,  for  example.Piper  (1974)  ).  Balas's  strategy  discovers 
complete,  feasible  solutions  which  successively  improve  the  objective  function 
value.  Tc  obtain  a complete  solution,  a partial  solution,  Jt,  is  iteratively 
augmented  by  assigning  a binary  value  to  some  variable  not  already  in  Jt.  The 
procedure  starts  with  an  'optimal'  but  infeasible  solution  and  attempts  to  reduce 
infeasibility  at  each  iteration.  If  it  is  determined  that  no  completion  of  can  be 
feasible  or  optimal  (and  after  a complete  solution  is  discovered),  a new  partial 
solution  is  defined  by  a special  backtracking  procedure  and  the  search  for  an 
improved,  complete  solution  continues.  The  strategy  enumerates,  either  im- 
plicity  or  explicity,  all  feasible  combinations  of  the  variables  so  that  it  guaran- 
tees an  optimal  solution. 

The  algorithm  developed  during  this  study  specializes  Balas's  strategy  to 

promote  efficiency  in  five  significant  ways: 

1.  Assigning  one  k-tour  to  level  1 in  Jt  requires  all  other  k-tours 
as  well  as  k-tours  to  be  assigned  to  level  0 to  invoke  equations 
[40]  and  [41].  Algorithm  logic  imposes  these  restrictions 
automatically,  no  feasibility  checks  (other  than  vehicle  network 
constraints,  if  any)  need  be  made  so  the  algorithm  essentially 
searches  only  for  improved  solutions. 
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2.  The  complementary  vehicle  relationships  in  equation  [41]  pro- 
mote efficiency  in  a manner  similar  to  the  cross  branching 
procedure  described  by  Martin  (1976). 

3.  Z is  treated  explicitly  as  a continuous  variable  so  it  is  not 
necessary  to  expand  the  problem  to  express  Z as  a function  of 
integer  variables. 

4.  The  single  variable  Z is  used  in  lieu  of  the  set  of  surplus  vari- 
ables employed  by  the  Balas  routine,  reducing  algorithm 
complexity. 

5.  Vehicle  and  vehicle-tour  set  memberships  are  defined  with  the 
objective  of  improving  rim  time  efficiency. 

A detailed  statement  of  the  algorithm  follows  the  discussion  below  which  pro- 
vides an  overview  of  the  logic  involved. 

Before  beginning  to  solve  Problem  01  at  iteration  ~ , the  'current  best 
solution'  to  01,  Z£,  is  known  to  be  UBK  as  defined  in  equation  ( 45 f and  only 
solutions  which  improve  Zc  need  be  examined.  The  algorithm  uses  an  efficient 
bookkeeping  procedure  to  invoke  equations  [40]  and  [41]  and  determine  which 
variables  might  best  augment  »L  at  the  j-iteration  within  01.  Before  iteration 
number  one,  all  tours  are  assigned  to  the  set  E^,  and  all  vehicles  to  set  Vc\ 
indicating  that  no  tours  have  arcs  associated  with  any  negative  y values  from 
equation  [37].  Each  tour  in  set  Et  has  the  same  coefficient,  in  all 

Benders'  cuts,  so  the  minimum  cost  T ^ is  the  best  k-tour  among  all  k-tours 
in  set  Efc.  This  one  k-tour,  the  'representative',  may  be  easily  identified  ,r 
is  the  only  one  from  E which  may  enter  solution  at  iteration”.  Any  other 
Tr|.  C E(_  would  necessarily  increase  Z over  the  lowest  attainable  value  u.-i  . 
the  representative  tour.  V1  is  defined  to  reduce  tour- membership  choc  inr  mo. 
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At  each  successive  ^-iteration,  Trk  Q Et  , which  employ  any  arc  associated 
with  a negative  y defined  by  LP  are  upgraded  from  Ej.  to  , the  set  of  tours 
which  are  candidates  for  use  in  01  during  the  x-iteration.  Representative  tours 
are  redefined  if  necessary. 

At  each  j-iteration  (within  01)  tours  eligible  to  augment  the  current  par- 
tial solution,  Jj.,  are  selected  from  the  list  of  candidates,  Ntc,  by  culling  any 
tours  which  can  be  shown  to  be  infeasible  or  non-optimal  by  Exclusion  Test  IL 
It  is  assumed  that  K'  vehicles  (K'  — K,  depending  upon  the  number  of  comple- 
mentarity constraints)  must  each  be  assigned  a tour  at  level  1 to  provide  a com- 
plete, feasible  network  design.  Tours  which  would  violate  network  constraints 
may  be  excluded  (assigned  to  set  Dt)  immediately.  Defining  WJ^.  as 


Wk  * l wttoc  Trk  tor  Trk  = Jt  • 

the  right-hand  side  of  each  equation  [43],  which  results  from  the  current  partial 
solution  is  B,,  + W'  . For  each  Benders'  cut  x,  the  minimum  cost  addition  to 
complete  a network  design  may  be  determined  knowing  T^  C Jt,  K ',  and 
vehicle  complementarity  relationships.  Each  tour  Trk  Q N.c  is  tested  in  each 
cut,  ac  , using  equation  [46]  to  determine  if  its  inclusion  [plus  the  minimum 
(Wrkx*  Wrk/c  ^ for  k ) required  to  complete  a network  design]  would  yield  a 
Z value  greater  than  the  current  best,  7 ! . Tours  which  could  not  improve  ZQ 
even  in  this  most  favorable  condition  are  excluded  by  assigning  them  to  set  Dt 
for  this  j-iteration.  A second  component  of  Exclusion  Test  ri  makes  a similar 
check  to  determine  if  J^,  Trk,  and  the  associated  minimum  cost  network  completion 
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would  exceed  the  cost  o!  using  the  commercial  air  mode  exclusively,  Bp,  de- 
fined by  Problem  LP.  Tours  which  could  not  yield  a network  cost  less  than  Bp 
in  this  most  favorable  condition  are  also  assigned  to  set  for  this  j-iteration, 
since  they  could  not  lead  to  a minimum  (system)  cost  solution  in  general.  In 
some  applications,  checking  time  may  be  reduced  by  identifying  the  'best'  tour 
of  a given  length  (number  of  arcs)  and  checking  conditions  of  step  l/B/2  and 
l/B/3  for  this  tour  first.  If  the  best  tour  of  the  given  length  fails  the  test,  all 
tours  of  the  given  length  may  be  excluded  immediately.  Tours  which  pass  the 
exclusion  test  are  eligible  to  augment  Jt  at  level  1. 

The  variable  to  enter  solution  is  determined  in  step  2 by  a heuristic  pro- 
cedure analogous  to  that  used  in  the  Balas  routine.  The  parameter  q forms  a 
linear  combination  (for  each  eligible  T ^)  of  the  sum  of  T^.  coefficients  in  all 
Benders'  cuts  and  its  maximum  coefficient  over  all  cuts.  Computational  exper- 
ience should  indicate  values  for  a which  provide  best  run  times. 

If  (K ' - 1)  vehicles  are  already  assigned  to  level  1 in  J^,  Exclusion  Test 
II  and  step  2/A  are  not  used.  Rather,  step  2/B  is  used  to  determine  the  feasible 
k-tour  which  minimizes  the  maximum  right  hand  side  of  all  Benders'  cuts,  yield- 
ing a Z value  of  Zp.  If  Zp  < Z£,  an  improved,  complete  solution  has  been  dis- 
covered and  is  recorded.  If  an  improved  network  cannot  be  defined,  the  back- 
tracking procedure  is  used  to  define  a new 

Jj  is  augmented  in  step  3 by  adding  the  entering  variable,  T p,  as  the 
rightmost  element  in  the  vector  as  is  done  in  the  Balas  procedure.  It  is 
(implicitly)  understood  that  all  other  candidate  T ^ (and  Tr]r  ) lor  k (and  k \ are 
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also  augmented  to  at  level  0,  between  the  old  elements  and  the  new  addition. 
This  is  accomplished  by  updating  the  set  membership  of  the  vehicles  involved  and 
and,  later,  by  the  logic  of  the  backtracking  step. 

The  backtracking  procedure  in  step  4 is  invoked  after  a complete  solu- 
tion has  been  discovered  or  when  it  is  determined  that  no  completion  of  Jt  can 
improve  Zc.  As  in  the  Balas  routine,  backtracking  sets  the  rightmost  positive 
element  of  Jt  to  0 (by  multiplying  the  element  by  -1)  and  dropping  negative  ele- 
ments (if  any)  to  the  right  of  this  element.  Set  membership  of  k (and  k ) is  up- 
dated appropriately  to  effect  changes  in  the  implicit  elements  of  Jt  at  level  0. 

The  current  best  solution  may  be  judged  optimal  in  the  backtracking  pro- 
cedure when  all  k-tours  (for  any  k)  are  included  as  0 explicitly.1  This  condition 
indicates  that  all  possible  combinations  of  including  all  k-tours  have  been  ex- 
amined (implicitly  or  explicitly)  and  further  reductions  in  Zc  could  be  attained 
only  by  excluding  vehicle  k,  an  infeasible  alternative.  If  this  k has  a comple- 
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mentary  vehicle,  k must  be  reassigned  to  the  set  of  candidate  vehicles,  Vt  , 

(k  is  assigned  to  a set  V of  vehicles  which  cannot  be  candidates  when  a k-tour 
enters  the  solution  at  level  1)  and  the  search  continued  until  both  k and  k (for 
some  k and  k ) have  all  candidate  tours  explicitly  included  in  Jt  at  level  0,  in- 
dicating that  all  possible  networks  using  either  k or  k have  been  (implicitly  or 
explicitly)  examined. 

1 Specifically,  if  no  other  vehicle  has  a tour  at  level  one  to  the  left 
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in  Jt. 


Step  0:  Initialization 


L 


A.  Initialize  iteration  k.  by  defining  the  'current  best  solution'  as 

Z0  - ub" 

and  set  Jt  = 0,  Vt  = 0,  V = 0,  V = 0,  NtC  = (N  - Et), 
v6  = i k | all  Trk  CEt|,  VtC=  (V  - V6),  V6  - 0,  j = 0. 

Complete  vehicle  networks  which  improve  ZQ  are  now  sought. 

B.  For  the  set  of  y — <0  defined  by  LP  at  iteration  ( k - 1),  some  of  the 

q/c 

arcs,  q,  may  not  have  been  associated  with  a negative  y 

value  from  a previous  iteration.  Remove  all  k tours  from  the  set  Et 

which  employ  any  of  the  q arcs  with  y — <0,  qCT  (k),  and  assign 

q*  Q 

them  to  the  set  of  candidate  tours,  NTtc.  If  k had  belonged  to  the  set  V6, 
remove  k from  Ve  and  assign  it  to  VC.  Go  to  step  1. 

Step  1;  Tour  Exclusion 
A.  Exclusion  I 

1.  If  Ej.  = 0,  go  to  step  1/B,  Exclusion  IL 

2.  For  each  vehicle  k which  has  tours  in  E , determine  the  tour  in  E. 

I » 

with  the  min  (Crk).  This  is  the  best  possible  k tour  in  E^  and  is 

called  the  representative  k-tour.  If  all  N tours  C Ej.,  STOP,  since 

the  vehicle  network  defined  by  the  representative  tours  is  optimal 

(but  not  practical  if  its  cost  is  greater  than  BQ).  Otherwise,  place 

each  representative  k-tour  in  set  Nt  . If  k C V°,  place  k in  set  \’£ 

since  the  representative  is  now  a candidate.  (If  the  null  k-tour  is 

defined,  it  will  always  be  the  representative  k-tour. ) 
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3.  Set  j = 1 and  go  to  step  2 since  «L.  = 0. 


B.  Exclusion  II 

1.  Increment  j (j  = j + 1).  If  only  one  more  vehicle  tour  is  required 
to  complete  a vehicle  network,  go  to  step  2/B. 

2.  For  each  k C vtc,  check  all  candidate  tours  (Trk  Q Nt°  ) to  deter- 
mine if  T , would  violate  a network  constraint  or  if,  for  cut 

rk 


K.  = 1,  2,  ...  (C 


1 T min  I W . . , W r,  l . [46] 

i'CVtc  LTrk'^NtC  rk>C  Fk  . 


W , >ZC-B  - W'  - 
rk k c K t/c  k 


k'?*k  Tr]r  C Nt°  (if  any) 

If  so,  remove  tour  rk  from  and  assign  it  to  the  set  of  excluded 
tours,  D^,  since  including  it  in  solution  would  either  be  infeasible  or 
would  increase  Z over  the  current  best  value,  ZQ.  If  this  is  the  last 
k-tour  C Ntc,  remove  k from  VtC  and  assign  it  to  V^. 

3.  For  each  k C check  all  candidate  tours  (T^  C N^)  to  determine 


C , >B„ 
rk  0 


* CrkTrk 
Trk  — ^t 


I 

min 

p 

{Crk"  Crk' 

k'CVtc 

LTrk’£Nt 

k'/k 

Trk' — Nt° 

(if  any) 

If  so,  remove  tour  rk  from  Ntc  and  assign  it  to  the  set  of  excluded 
tours,  Dt>  since  including  it  would  increase  Z over  BQ.  If  this  is  the 
last  k-tour  C Nj.c,  remove  k from  \'£  and  assign  it  to  Vcl. 
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4.  The  set  of  tours  T , C N.  for  vehicles  kC  V,  are  eligible  to  aug- 

rk  i ~~  1 

ment  Jt  at  level  one.  If  Ntc  = 0,  no  completion  of  Jt  will  improve  Zc; 
set  j = j + 1 and  backtrack  by  going  to  step  4.  Otherwise  go  to  step  2. 


Step  2:  Determine  variable  to  enter  at  level  one.  If  only  one  tour  is  required 
to  complete  the  vehicle  network,  go  to  step  2/B;  otherwise,  step  2/A. 

A.  Entering  tour  will  not  complete  the  vehicle  network. 

1.  The  tour  to  enter  solution  at  this  iteration  is  the  tour,  T *,  , for  some 

rk 

c c 

k C V and  some  T , C N.  which  provides  the 
t rk  t 

+ (1  - a ) 


min  ( a 
Trk=N<c 


K. 

1 XV 

K.  — 1 


vkic 


max  W , 

t rkx 


in  which  a is  a parameter  0 — a— 1.  If  a - 0,  the  tour  which  adds 

the  least  to  any  cut  would  enter;  if  a - 1,  the  tour  which  has  the 

smallest  sum  of  cut  coefficients  would  enter. 

2.  Go  to  step  3 to  augment  with  the  tour  T^  . 

B.  One  tour  is  required  to  complete  the  vehicle  network. 

1.  For  eligible  vehicles  (k  C V,C  ) select  the  tour  T * from  the 

t rk 

Q 

T , C N+  which  violates  no  network  constraints  and  which  pro- 

rk  _ t 

vides  the 


min  Z = min  1 max 

p Th.cn,'  I ' 

2.  If  Z > Z or  if  no  T , is  feasible  with  respect  to  the  network  con- 
p c rk 

straints,  go  to  step  4 to  backtrack  since  no  completion  of  can 
improve  the  value  Z . 


(B* 


W’  ) 
tx 


+ W 


rk'cj  ( 
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3.  If  Z ^ , go  to  step  3 to  augment  with  the  tour  T 


Step  3:  Augment  Jt  by  adding  T^k  at  level  one. 

A.  Updates 

1.  Assign  as  the  rightmost  element  of  Jt. 

2.  Update  W ' V/c  and  C ' . 

tic  t 

3.  Remove  k from  Vtc  and  assign  it  to  V^.  If  k C Vtc  or  if  k C V0*, 

■ — Q 

remove  k and  assign  it  to  V . Remove  T*  from  Nt  . 

t rk  L 

4.  Assign  all  k C to  Vc  and  set  = 0. 

Assign  all  T C Dt  to  and  set  = 0. 

B.  Continue  implicit  enumeration. 

1.  If  the  vehicle  network  is  not  complete,  go  to  step  1/B. 


2.  If  the  vehicle  network  is  complete,  Z will  be  — Z . Record  the 

P 


solution  and  set  Z„  = Z . Go  to  step  4 to  backtrack. 

c p 


Step  4;  Backtrack 

R 

A.  Find  the  rightmost  positive  tour  in  the  vector  Jt,  • If  there  are 
negative  tours  to  its  right,  delete  them  from  J.  and  reassign  them  to 
Ntc.  If  the  associated  k had  belonged  to  V^,  reassign  it  to  Vtc. 

B.  If  there  are  candidate  k-tours  not  in  Jj.,  go  to  step  4/C.  If  all 

candidate  k-tours  are  in  Jj.  and  if  there  are  no  tours  at  level  one  to 

the  left  of  T^,  go  to  step  4/D.  If  all  candidate  k-tours  are  in  Jt  and 

R R 

if  there  is  some  tour  at  level  one  to  the  left  of  Trk,  multiply  Trkby  -1 
and  go  to  step  4 /A. 
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C.  Backtrack,  Make  T , negative  to  specify  that  the  tour  is  included  in 

rk 

Jt  at  level  zero.  Remove  k from  Vt  and  assign  it  to  V^.  If  T;  C V^.  , 
remove  it  and  assign  it  to  Vtc.  Go  to  step  1/B  to  begin  another  iteration. 

D.  Optimality  Conditions.  If  k has  no  complementary  vehicle,  STOP;  since 
the  optimal  solution  has  been  identified  (all  k tours  have  been  fathomed). 
If  k C Vt*  also;  STOP,  and,  again,  the  optimal  solution  has  been  identi- 
fied. If  k (£  Vtc,  assign  k to  V^,  remove  k from  V and  assign  it  to 
make  T ^ negative  to  specify  that  the  tour  is  included  in  at  level  zero, 
and  go  to  step  1/B  to  begin  another  iteration. 
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CHAPTER  V 


TOUR  CONSTRUCTION 

In  many  applications,  particularly  in  cases  involving  a large  number  of 
points,  it  may  be  difficult  to  predefine  a limited  number  of  tours  to  make  avail- 
able to  Problem  01.  In  fact,  in  most  applications,  it  may  be  preferable  to  con- 
struct tours  using  an  appropriate  model  formulation  and  solution  algorithm. 
Possibilities  for  tour  construction  are  examined  in  this  chapter. 

A basic  requirement  of  any  tour  construction  approach  is  the  definition 
of  some  network  over  which  vehicle  tours  may  be  determined.  Figures  9 and  10 
depict  characteristics  of  one  type  of  network  which  may  be  used  to  define  tours. 
A simple  case  involving  two  complementary  vehicles  and  a set  of  points  is  shown 
in  Figure  9.  Each  actual  point,  i,  appears  in  each  vehicle  subnetwork  coded  as 
point  ik.  Should  more  than  one  visit  be  allowed  to  point  i in  subnetwork  k,  a 
number  of  'dummy'  k-points  need  be  defined  to  represent  the  possibility.  Fig- 
ure 10  indicates  a device  used  in  the  formulation  to  constrain  the  number  of 
visits  to  each  ik  point  to  be  at  most  one.  In  general,  the  network  to  define  tours 
may  be  formulated  as  a shortest  path  problem,  or,  as  shown  in  Figure  9,  as  a 
flow  circulation  problem.  Since  the  number  of  vehicles  required  to  define  a 
network  is  one  (K'  = 1)  in  this  example,  one  unit  of  'flow'  circulates  through 
the  network  describing  either  a k or  a k-tour  from  dummy  origin,  O,  to  dummy 
termination,  T,  and  returning  to  ()  along  the  return  arc  (T,  6).  The  network 
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FLOW  CIRCULATION  NETWORK  FOR  VEHICLE  TOUR  CONSTRUCTION 


igure  9 


FLOW  RESTRICTION  THRU  POINT  i 


concept  may  allow  a variety  of  logical  network  constraints  to  be  defined  to  model 
actual  requirements  in  a particular  application. 

The  distinction  between  the  multiple- vehicle  routing  problem  and  the 
transportation  system  planning  problem  is  made  evident  in  the  flow-circulation  net- 
work. In  the  former  problem  type,  each  (actual)  i-point  is  represented  in  each 
subnetwork  by  a single  ik-point  and  each  i-point  must  be  visited  once  and  only 
once  in  all  K ' subnetworks  combined.  If  only  one  vehicle  be  used,  a traveling 
salesman  problem  must  be  solved  to  determine  the  minimum  cost  flow  (path) 
from  6 to  T which  is  constrained  to  visit  all  ik-points.  In  the  transportation 
system  planning  problem,  not  all  ik-points  need  be  visited  by  a feasible  k-tour, 
and  some  i-points  may  be  visited  more  than  once.  Hence,  the  minimum  cost 
flow  path  from  to  is  sought,  and  the  points  visited  on  the  tour  are  not 
predetermined. 

Three  strategies  by  which  the  tour  construction  approach  may  be  imple- 
mented are  examined  in  this  chapter.  An  implicit  enumeration  approach  is  de- 
scribed first.  An  alternate  formulation  permits  use  of  efficient  LP  procedures 
in  a branch-and-bound  strategy  discussed  in  the  second  section.  Finally  a strat- 
egy to  construct  tours  by  solving  a simple  flow  circulation  problem  is  considered. 


An  Implicit  Enumeration  Approach 

The  primary  (binary)  variable  in  any  tour  construction  approach  must  be 


one  defined  as  a : 


1 if  vehicle  k traverses  the  arc  from  point  i to  point  j 
(i,  j,  and  k are  implicitly  related  by  the  index  q), 

0 otherwise. 


i 4 
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This  variable  bears  a close  relationship  to  the  binary  variables  used  in 


Problem  PI  since 


[47] 


for  the  appropriate  k.  The  cost  of  a vehicle  tour  may  be  defined  in  detail  using 


the  new  variables: 


q£Tq(k) 


cd% 


Ck  I1  - Tok> 


in  which  the  costs  of  traveling  from  i to  j and  of  servicing  point  j are  included 
in  Cq  and  fixed,  investment-type  expenses  for  vehicle  k are  represented  by  C^. 

Combining  the  flow  circulation  concept  and  the  new  cost  function,  the 
transportation  system  planning  problem  may  be  formulated  as  shown  in  Figure  11. 
Problem  P2  is  equivalent  to  Problem  PI,  which  was  examined  in  the  previous 
chapter. 

Material  flow  constraints  given  in  equations  [49]  to  [52]  are  the  same  as 
those  used  in  Problem  PI  with  one  exception.  The  right  hand  side  of  each  vehicle 
arc  capacity  constraint  (equation  [51]  ) is  now  defined  with  respect  to  a rather 
than  Trl<  variables  using  the  relationship  noted  in  equation  [47]  above. 

The  remaining  constraints — equations  [53]  through  [59] — impose  the  flow 
circulation  concept  to  construct  tours.  Practical  limitations  on  tour  length  are 
represented  by  equation  [53].  In  some  applications,  more  complicated  constraints 
may  be  required,  for  example,  to  place  lower  limits  on  tour  length  as  well,  or  to 
allow  tour  length  to  fall  in  one  of  several  feasible  ranges.  Equations  [54]  and  [55] 
allow  each  ik-point  to  be  visited  at  most  once,  imposing  the  form  shown  in 


74 


Problem  P2: 


Min 


msd 


L 

2 C.  Po(t)  + 2 C rA  , S rA , + 2 C rT  . S rT  . 

£ = 1 * msd  msc*  msd  msd  msd  msd 


+ 


Q 

V* 

q~i 


K 

C a + 2 

* k=l 


CF  (1 


Tok> 


[48] 


Subject  to: 


Ship  Time  Standards 


2 Tlt,fmsdPt'+  2 {~~>CA  ^ 

fe|i(m)  msd 


oin  ,1  qTT  . 

£ ‘msd  ‘A  Amsd  ' 'msd  msd  Tmsd  "msd  m* 


m-1.  2. 


M 


Transshipment  Point  Capacity 
L 


£-1 


*Ufm.d  pl  s bl  • '“‘.2 ' 


[491 


[50] 


Vehicle/arc  Capacity 
L 


aq£  ‘msd 


*msd  ^q  aq  > *1  1, 2,  . . . , Q 


[51] 


Ship  Requirements 


("t ) r A r T 

8msd,i.  !>A  * Smsd  + Smsd  ~ 1,0  ’ 

' “1 


V msd 


[52] 


Tour  Time/Distance  Length  Limit 


2 d a s (]  ; 

q^Tq(k)  q q 


k-1,  2, 


1 53] 


Figure  11. --Alternate  Formulation  of  the  Transportation  System 

Planning  Model 
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Problem  P2:  (Continued) 


Flow  Into  Point  ki 


2 aQ  - Yik  = ° ; v ik 
qe  Ej(k)  q 


Flow  Out  of  Point  ki 


Y.,  - 2 a = 0 ; V ik 

lk  qeBj(k)  1 


Complementary  Vehicle  Relationships 


TQk  + Tjj;  = 1 ; for  appropriate  k,  k 


Vehicle  Network  Constraints 

As  required  in  a specific  application  to  define  the  most 
appropriate  network  flow  problem,  and  to  impose  the 
logical  structure  of  the  network;  i.  e. , 


- 0 ; k-1,  2, . . . , K 


q6B0(k)  q 


Lower/ Upper  Bounds  on  Flows 


0 — aq  — 1;  V q 


0 - Yilc  - 1;  V ik 
0 ^ Tok  ^ 1;  V k 
Non-negativity  and  Integer  Constraints 

'f  ' Smsd’  Smsd  2 0 ' «•  msd 

*q  • Yik-  Tok  »•  1 '•  T 1.  k 


Figure  11  (Continued) 


J 


Figure  10.  If  an  arc  incident  to  point  ik  is  placed  on  the  k-tour,  these  con- 
straints assure  that  some  arc  leaving  point  ik  will  also  be  placed  on  the  k-tour. 

A connected  path  from  Oj^  to  Tk  will  therefore  be  defined.  No  subtours  allow- 
ing flow  to  circulate  amongst  a set  of  k-points  not  on  this  path  will  be  permitted, 
since  that  additional  flow  would  only  increase  costs  and  would  therefore  not  be 
optimal.  (This  type  of  problem  is  encountered  by  the  subtour  elimination  ap- 
proach to  the  traveling  salesman  problem  and  is  discussed  in  Chapter  II.  ) 

Complementary  vehicle  relationships  may  also  be  used  in  this  formula- 
tion as  shown  by  equation  [56],  Additional  network  constraints  important  in  a 
particular  application  may  also  be  added  according  to  the  ingenuity  of  the  analyst. 
In  particular,  logic  required  to  use  the  flow  circulation  strategy  must  be  invoked 

by  constraints  such  as  equation  [57],  which  essentially  allows  the  variable  T . 

ok 

to  be  defined,  indicating  whether  or  not  the  null  k-tour  is  used  in  solution. 

Lower  and  upper  bounds  for  flows  on  each  arc  in  the  conceptual  network  are 
stated  in  equation  [58],  and  correspond  to  the  concepts  described  earlier. 

Problem  P2  is  also  a large-scale,  mixed  01  integer,  linear  program- 
ming problem  which  may  be  reduced  to  two  more  simple,  interacting  problems 
by  aPPtying  the  concepts  of  Benders'  decomposition.  The  resulting  linear  pro- 
gram, Problem  LP',  is  stated  in  Figure  12.  LP'  is  the  material  flow  component 
and  is  the  same  as  Problem  LP  with  the  exception  that  the  right  hand  sides  of  the 
vehicle  arc  capacity  constraints  (equation  |(i3|  ) are  now  stated  in  terms  of  the 
a(j  instead  of  the  Tr|<.  A particular  network  design  is  communicated  to  LP'  by 
the  set  of  values  assigned  to  the  binary  variables  a(j.  In  turn,  LP'  feeds  back 
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Problem  LP': 


Min  f 2 C P«^)  + 2 C rA  . S rA  , + 2 C S fT  1601 

msd  2=1  L £ msd  msd  msd  msd  msd  msd 
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Ship  Requirements 
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2 e pil>  + s +s 
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V msd 
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m 1 1 


/3i  (62] 


yq  [63] 


Smscll  [fi4l 


Non-negativity  Constraints 


p/>  . sit. . S 11 


\ > ®msd  > Smscj  - 0 ; V-£, , msd 


Figure  12. —Problem  LP’ 


information  concerning  the  efficacy  of  the  network  design  by  adding  a constraint 


to  Problem  MP  using  the  optimal  values  of  the  associated  dual  variables. 

The  interacting  integer  program,  Problem  MP,  which  consists  of  the 
flow  circulation  constraints  augmented  by  the  Benders'  cuts,  is  stated  in 
Figure  13.  As  in  Problem  01,  the  are  non-positive  and  (as  defined  by 
equation  [38])  is  non-negative.  Since  a feasible  solution  to  LP'  is  guaranteed  by 
the  assumption  that  exclusive  use  of  commercial  air  transportation  is  feasible, 
only  cuts  of  the  form  in  equation  [66]  are  necessary. 

Problem  MP  makes  evident  a network  structure  which  may  be  used  to 
develop  an  efficient,  implicit  enumeration  algorithm  based  on  the  Balas-type 
strategy.  As  in  the  algorithm  developed  to  solve  Problem  01,  the  theory  of 
Benders’  decomposition  provides  an  initial  upper  bound  on  the  solution  to  MP. 
Furthermore,  most  constraints  may  be  invoked  by  algorithm  logic  rather  than 
time  consuming  search;  the  variable  Z may  be  used  to  evaluate  the  optimality7 
of  a solution,  treating  it  as  a continuous  variable  and  circumventing  the  need  to 
add  a surplus  variable  to  each  Benders'  cut.  Tour  length  constraints  may  be  in- 
voked easily.  Branching  may  be  accomplished  according  to  the  set  of  vehicle 
arcs  which  depart  from  the  last  point  on  a partial  tour,  which  may  be  constructed 
by  beginning  at  point  O in  Figure  9 and  progressing  toward  point  f in  each  vehicle 
subnetwork. 
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Problem  MP 


Min  Z 
Subject  to: 

Benders’  Cuts 


Q 


K 


Z ^ 2 ( C - bn  y /r  ) a + 2 C (1  - T ) + B ; 

q=1  ' q q 1 qK  1 q k=1  k ok  K 


K.  =1,  2,  . . . , K 


Tour  Time  Distance  Length  Limit 


d a — H,  ; k-1,  Z, . . . , K 

q q k 


q£Tq(k) 

Flow  Into  Point  ki 


2 a - Yik  = 0 ; V ik 

qeEjjk) 

Flow  Out  of  Point  ki 


Y. 


ik 


2 a = 0 ; V ik 

q 


qe  Bj  (k) 

Complementary  Vehicle  Relationships 

Tok  + ^ok  = 1 ’ f°r  aPProPr^ate  k 

Vehicle  Network  Constraints 

As  required  in  a specific  application  to  define  the  most 
appropriate  network  flow  problem,  and  to  impose  the 
logical  structure  of  the  network;  i.  e. , 


rok“  2 aq  0 5 k 1,2 K 

q£B0(k) 


Figure  13. — Problem  MP 
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A Branch-and-Bound  Approach 

The  column  generation  technique  developed  in  Chapter  III  to  solve  the 
material  flow  problem  is,  from  one  viewpoint,  merely  one  way  to  define  the 
elements  in  the  column  associated  with  the  next  (continuous)  variable  to  enter 
solution  in  the  linear  program.  Another  view  of  the  procedure  reveals  that  it  is 
an  efficient  way  of  assigning  values  to  a large  number  of  binary  variables:  a 

qt 

tj^,  and  X.^.  Since  this  set  of  binary  values  must  satisfy  a set  of  feasibility 
conditions,  the  column  generation  procedure  solves  an  integer  program  in  the 
course  of  determining  the  level  at  which  a continuous  variable  will  enter  solution 
in  a linear  program.  With  this  idea  in  mind,  it  seems  that  an  efficient  tour  gen- 
eration procedure  might  be  developed  by  formulating  an  appropriate  model  and 
Using  simple  subproblems  (such  as  shortest  path  problems)  to  generate  tours 
without  solving  a large  integer  program  directly.  One  model  formulated  in  the 
attempt  to  effect  this  strategy  appears  in  matrix  form  in  Figure  14. 

Constraints  of  the  type  in  equation  [75]  impose  vehicle  network  logic-type 
constraints,  complementary  vehicle  relationships,  and  tour  length  limits.  Equa- 
tion  [76]  constrains  each  tour  to  visit  each  point  at  most  once  (multiple  visits  could  be 
permitted  by  defining  another  'dummy'  point  to  represent,  for  example,  the 
second  visit  to  one  point).  The  entire  model  represents  the  integer  subproblem 
resulting  from  Benders'  decomposition  and  equation  [77]  imposes  the  Benders' 
cuts.  Finally,  equation  [78]  states  that  each  vehicle  may.  traverse  at  most  one 
tour.  Relaxing  the  requirement  that  the  T^  need  by  binary,  the  problem  becomes 
a linear  program  which  has  exactly  the  structure  of  the  material  flow  problem. 
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Equations  [78]  are  of  the  Generalized  Upper  Bounding  type,  logic  con- 
straints are  analogous  to  the  ship  time  standards,  equations  [76]  are  essenti- 
ally the  same  as  the  vehicle  arc  capacities,  and  Benders'  cuts  correspond  to 
transshipment  point  capacity  limitations.  The  GUB  and  column  generation  pro- 
cedures described  in  Chapter  III  could  be  applied  directly  to  this  formulation. 

Tours  would  be  generated  as  desired;  i.  e. , as  a by  product  of  determining  the 
entering  variable.  The  difficulty  posed  by  this  formulation  is  that  the  program 
could  assign  some  level  to  Tr^  between  0 and  1 so  that  several  k-tours  might  be 
in  the  optimal  solution  of  the  linear  programming  problem  simultaneously,  each 
at  some  fractional  value.  Such  a solution  would  have  no  practical  meaning. 

However,  this  model  and  approach  would  be  efficient  timewise.  The 
relaxed  version  of  integer  problems  is  commonly  used  to  give  some  indication 
of  the  solution;  or,  at  least,  to  provide  a bound  for  the  optimum  value  of  the  ob- 
jective. Binary  values  would  not  be  guaranteed  to  result  from  the  linear  pro- 
gram since  the  constraint  matrix  is  not  unimodular  (unimodularity  is  shown  to 
be  a sufficient  condition  for  integer  solutions  to  result  from  application  of  the 
linear  programming  procedure  by  Hu  (1969)  ).  However,  computational  exper- 
ience might  show  that  integer  solutions  are  obtained  routinely.  Alternately, 
solutions  which  might  be  rounded  (say,  by  some  heuristic  procedure)  to  'very 
good'  integer  solutions  might  be  attained. 

In  fact,  it  is  possible  to  improve  computational  efficiency  even  further 
by  application  of  an  additional  large-scale  programming  technique.  Blocks  of 
constraints  each  representing  limitations  on  a specific  vehicle  may  be  identified 
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in  the  model  presented  in  Figure  14.  The  set  of  Benders'  cuts  and  certain 


logical  relationships  which  relate  several  vehicles  may  be  viewed  as  coordin- 
ating constraints  relating  vehicle  tours.  This  structure  may  be  solved  effici- 
ently by  Dantzig-Wolfe  decomposition,  or  by  an  extended  version  of  the  GUB 
method  which  treats  such  block-diagonal  forms  (see  Lasdon,  1970).  Both  of 
these  methods  would  maintain  a B-^  for  each  vehicle  block,  rather  than  for 
the  entire  constraint  set  and  would  thereby  improve  storage  and  run  time  effic- 
iency. At  each  iteration,  a new  tour  for  a single  vehicle,  rather  than  a set  of 
tours  for  all  K'  vehicles,  would  be  constructed  by  the  column  generating  pro- 
cedure. In  any  event,  the  approach  might  be  used  to  generate  tours  used  in  the 
first  few  /c-iterations  to  develop  a set  of  Benders'  cuts  which,  upon  application  of 
the  integer  programming  algorithm,  might  lead  to  an  optimal  solution  rather  quickly. 

Alternately,  and  preferrably,  the  linear  integer  program  in  Figure  14 
could  be  solved  using  a branch-and-bound  procedure  such  as  that  devised  by  Land 
and  Doig  (1960),  and  Dakin  (1965).  Forrest,  Hirst,  and  Tomlin  (1974)  describe 
efficient  techniques  for  solving  large-scale  mixed,  integer  linear  problems  of 
this  type.  The  branch-and-bound  strategy,  in  combination  with  the  efficient  pro- 
cedures applicable  to  the  linear  program,  appears  to  offer  an  efficient  approach 
for  solving  this  integer  problem  directly. 

Tour  Construction  Via  Flow  Circulation 

Problem  MP  provides  an  (apparently)  compact  model  for  solution  by 
implicit  enumeration  since  it  identifies  feasibility  conditions  which  might  be 
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checked  rapidly.  In  addition,  the  model  suggests  the  opportunity  to  construct 
tours  using  very  efficient  network  flow  subproblems. 

The  set  of  feasible  solutions  to  Problem  MP  (and  hence  the  optimal 
solution)  is  a subset  of  the  feasible  solutions  to  the  straightforward  network 

flow  problem  (i.  e. , ( an  T )C  A ) defined  by  equations  [68],  [69],  and  [72]. 

— 4 -ok  

A procedure  which  searches  efficiently  among  these  extreme  points  would  there- 
fore appear  to  be  well  worth  evaluating.  Several  possibilities  exist. 

A modest  amount  of  research  which  develops  procedures  for  invoking 
constraints  by  formulating  them  as  components  of  the  objective  functions  has 
been  reported  in  the  literature:  Anthonisse  (1973),  Bradley  (1971),  Glover  and 
Woolsey  (1972),  Kendall  and  Zoints  (1977).  The  procedures  all  remove  an  equa- 
tion from  the  constraint  set,  multiply  it  by  some  (numerically  large)  integer  and 
incorporate  the  resulting  function  into  the  objective.  Some  of  these  procedures 
operate  on  constraints  two  at  a time,  a policy  which  could  be  implemented 
effectively  for  the  current  problem,  since  Benders'  cuts  are  added  one  at  a 
time.  The  integer  multipliers  must  satisfy  certain  conditions  (particularly  that 
they  be  relatively  prime  with  respect  to  each  other)  to  assure  that  the  under- 
lying constraints  are,  in  fact,  satisfied  by  the  optimal  solution  to  the  substituted 
problem.  The  procedures  are  distinguished  by  the  conditions  each  establishes 
for  selecting  specific  values  for  the  integer  multipliers.  The  difficulty  which 
none  of  these  procedures  has  overcome  is  that  the  multipliers  must  increase  in 
numerical  value  for  each  successive  constraint  incorporated  into  the  objective. 
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Consequently,  the  multipliers  must  be  assigned  enormous — in  fact  computation- 
ally prohibitive— values  when  more  than  a few  constraints  are  involved. 

This  approach  would  be  useful  if  the  optimal  solution  to  Problem  P2 
could  be  identified  after  augmenting  only  a few  Benders'  cuts  to  MP  (in  the 
absence  of  tour  constraints  such  as  equations  [53]  ).  It  is  expected  that  the 
number  of  cuts  required  to  identify  the  optimum  in  problems  of  modest  size 
would  be  on  the  order  of  10-15.  This  is  still  more  than  the  four-to-six  con- 
straints which  may  be  handled  effectively  by  available  procedures  of  this  type. 
However,  this  tactic  would  be  successful  for  solving  MP  at  each  of  the  first 
few  k. -iterations.  After  relegating  constraints  to  the  objective  function. 

Problem  MP  reduces  to  a simple  network  flow  problem  for  which  a number  of 
very  efficient  algorithms  are  available  (for  example,  see  Fulkerson  (1961), 
Durbin  and  Kroenke  (1967),  and  Bradley,  Brown,  and  Graves  (1976)  ). 

An  alternate  approach  would  apply  the  Lagrange  multiplier  technique  to 
incorporate  equations  [66]  and  [67]  into  the  objective  function.  Everett  (1963) 
published  the  first  application  of  Lagrange  multipliers  in  solving  integer  pro- 
grams. Brooks  and  Geoffrion  (1966)  developed  an  efficient  means  of  imple- 
menting Everrett's  approach  using  linear  programming.  Nemhauser  and 
Widhelm  (1971)  formalized  and  extended  this  work.  Others,  including  Geoffrion 
(1974)  and  Shapiro  (1971), have  contributed  to  the  development  of  this  approach 
to  integer  programming. 

Application  of  the  Brooks/Geoffrion  approach  to  Problem  MP  yields 
Problem  BG  in  Figure  15.  The  objective  of  BG  is  to  minimize  a convex  com- 
bination (see  equation  1 791  ) of  the  Z (objective  values  in  MP)  which  result  from 
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Problem  BG: 


J 

Min  Z (ZJ)  A..  [79] 

i — 1 


Figure  15.  — Problem  BG 


88 


/ 


each  of  the  J feasible  extreme  points  (vehicle  networks)  of  MP;  that  is  for 
each  ( a^,  T^)  C A,  where  the  set  A is  defined  by  equations  [68]  through  [73], 
Equations  [81]  and  [82]  assure  that  the  A.  multipliers  generate  solutions  which 
are  feasible  with  respect  to  equations  [66]  and  [67]  in  MP.  BG  is  a linear  pro- 
gram in  the  continuous  variables.  The  purpose  of  BG  is  to  generate  optimal 
values  of  the  Lagrange  multipliers  for  MP.  These  multipliers  are,  in  fact,  the 
set  of  pQ,  , and  "p^  variables  which  are  the  Simplex  multipliers  in  BG. 

The  relationship  of  this  approach  to  the  current  line  of  investigation  is 
evident  upon  considering  the  criterion  which  determines  the  entering  column 
at  each  iteration  of  the  Simplex  method  applied  to  BG.  The  criterion  to  identify 
the  best  column  to  enter  solution  (min  ”C.  - Z."<  0)  is 

J J J 
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which  may  be  expressed  for  closer  scrutiny  as 
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This  criterion  evaluates  the  j**1  extreme  point,  (2q’  lok)-A  • ^ a 
candidate  to  enter  solution  to  the  linear  program,  BG,  given  the  current  values 
of  the  Simplex  multipliers:  pQ,  p , andp^.  The  constant  term,  p ^,  and  all 
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terms  involving  and  may  be  viewed  as  defining  the  objective  function 
for  a network  flow  problem  defined  on  the  set  of  extreme  points,  A . Over- 
looking one  difficulty  which  is  discussed  below,  this  approach  would  solve  MP 
very  efficiently  by  solving  BG,  a linear  program,  using  a network  flow  algorithm 
to  generate  entering  columns  as  in  equation  [84].  At  the  optimal  solution  to  BG, 
optimal  values  of  the  Lagrange  multipliers  ( pQ,  p^,  and  pk  ) are  known.  The 
vehicle  network  generated  by  these  multipliers  is  the  optimal  solution  to  MP 
at  this  k.  -iteration. 

The  difficulty  which  hampers  this  tactic  is  the  term  involving  z)  in  equa- 
tion [84].  As  noted  in  equation  [66],  is  determined  as  a function  of  the 
( a j , T u ) vector  evaluated  by  all  k Benders'  cuts.  It  is  therefore  not  a term 
which  may  be  related  simply  to  a single  arc  in  the  network  flow  problem  as  all 
other  terms  in  equation  [84]  may  be.  There  is,  therefore,  apparently  no  simple 
network  flow  subproblem  which  may  be  used  in  the  required  column  generation 
procedure.  The  summation  involving  the  p^  could  be  restricted  to  equal  1.  0 and 
cause  the  ZJ  term  to  be  eliminated  so  that  the  network  flow  subproblem  could  be 
used  as  a column  generator.  However,  research  efforts  have  not  identified  con- 
ditions which  support  this  as  a legitimate  restriction. 

A second  difficulty  may  be  encountered  using  the  Lagrange  multiplier 
approach.  Problem  BG  may  have  an  optimal  solution  and  therefore  a set  of 
associated,  'optimal'  I^agrange  multipliers.  But  these  multipliers  might  not 
generate  the  optimal  vehicle  network  for  Problem  MP.  This  'gap'  phenomenon 
is  incurred  when  it  is  not  possible  for  the  Lagrange  multipliers  to  generate  the 


90 


right  hand  sides  (resources)  of  the  associated  constraints.  In  such  cases,  the 
approach  fails  to  define  the  optimal  solution  to  MP.  The  literature  does  not 
record  computational  experience  which  would  allow  one  to  develop  an  expecta- 
tion of  the  likelihood  of  the  occurrence  of  'gaps'  in  the  current  problem. 
Nemhauser  and  Widhelm  (1971)  discuss  this  problem  at  length  and  suggest  a 
means  of  detecting  gaps.  Their  approach  is  more  sophisticated  than  the 
Brooks/Geoffrion  method  applied  above,  but  it  introduces  nonlinear  relation- 
ships which  would  complicate  solution  of  the  transportation  system  planning 
problem. 

In  any  event,  BG  could  be  used  to  solve  MP  if  the  criterion  in  equation 
[84]  could  be  defined  by  an  efficient  procedure.  In  the  absence  of  an  optimizing 
procedure,  the  formulation  could  be  the  basis  of  a number  of  heuristic  proce- 
dures which  could  be  devised  to  define  an  objective  for  the  network  flow  sub- 
problem. Since  problems  of  this  type  may  be  solved  very  efficiently,  even 
heuristics  which  require  solution  of  a large  number  of  network  flow  problems 
(in  identifying  the  best  MP  solution  at  the  k -iteration)  would  be  practical  and, 
perhaps  preferable  to  the  implicit  enumeration  procedure  described  in  the  pre- 
vious subsection.  A heuristic  procedure  which  identifies  'good'  solutions  may 
be  useful  in  generating  'good'  Benders'  cuts  using  LP'.  That  is,  even  though 
the  network  supplied  to  LP'  is  not  truly  optimal  in  MP  at  the  current  k. -iteration, 
it  would  allow  LP'  to  furnish  additional  information  to  MP  (in  the  form  of  the 
Benders'  cut  added  at  the  next  iteration)  to  progress  toward  the  optimal  solution. 
The  optimizing  routine  using  the  implicit  enumeration  strategy  may  then  be 
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required  for  the  last  few  iterations.  Total  run  time  may  be  reduced  appreciably, 
but  computational  experience  is  necessary  to  estimate  computation  time  required 


for  problems  of  realistic  size. 


T 
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CHAPTER  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 
FOR  FUTURE  RESEARCH 

Conclusions 

The  transportation  system  planning  models  developed  during  this  study 
appear  to  be  capable  of  incorporating  a number  of  considerations  important  to 
the  design  of  actual  systems.  In  general,  this  study  contributes  to  early  work 
on  the  class  of  problems  which  require  a network  to  be  synthesized  to  satisfy,  in 
an  optimal  manner,  a set  of  conditions  concerning  the  network  operation.  A 
large  number  of  scheduling  and  vehicle  routing  type  problems  may  be  shown  to 
be  members  of  this  class. 

In  particular,  the  formulations  presented  in  this  report  are  capable  of 
imposing  important  operating  requirements  concerning  ship  time  performance 
which  may  be  used  to  evaluate  the  potential  effectiveness  of  a transportation 
plan  in  a large-scale  logistics  network.  In  addition,  the  approaches  developed 
are  capable  of  solving  problems  which  involve  symmetric  distance  matrices,  a 
condition  which  precludes  the  use  of  most  procedures  designed  to  solve  the  more 
simple,  multiple-vehicle  routing  problem. 

It  is  expected  that  the  material  flow  problem  described  in  Chapter  III 
will  be  useful  in  a number  of  practical  applications,  even  though  it  does  not 
provide  a mathematically  optimal  system  design.  This  expectation  is  based  on 
the  observation  that  much  of  the  recent  work  in  operations  research  is  directed 
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toward  providing  tools  of  analysis  which  may  be  used  to  augment  human  judgment 
rather  than  replace  it  with  inflexible  requirements  imposed  by  mathematical 
necessity.  In  complex  logistics  systems,  it  may  be  difficult,  if  not  impossible, 
to  make  explicit  all  constraints  and  operating  considerations  important  to  the 
ultimate  system  design.  Use  of  the  material  flow  model  permits  this  type  of 
flexibility  while  providing  an  efficient  solution  procedure  based  on  concepts 
known  to  be  efficient  computationally  in  prior  studies. 

The  optimizing  method  presented  in  Chapter  IV  provides  the  capability 
to  determine  a mathematically  optimal  solution  in  problems  of  modest  size,  and 
may  be  used  to  indicate  a reference  system  design  which  might  be  further  modi- 
fied as  desired  by  managerial  judgment.  The  approach  does  not  appear  to  be 
severely  restricted  by  the  need  to  identify  likely  tours  for  inclusion  in  system 
design,  since  it  is  expected  that  such  tours  could  be  readily  identified  in  most 
practical  cases.  The  algorithm  which  enumerates  all  feasible  tours  extends 
capability  to  solve  all  problems  of  modest  size.  The  proposed  solution  approach 
which  includes  use  of  Benders'  decomposition  and  an  implicit  enumeration  algo- 
rithm based  on  the  Balasian  approach  is  expected  to  be  efficient,  since  both  con- 
cepts have  proven  computationally  attractive  in  numerous  prior  applications. 

The  greatest  potential  for  algorithm  autonomy  is  presented  by  the  tour 
construction  approaches  described  in  Chapter  V.  Identification  of  the  material 
flow-circulation  formulation  of  the  network  design  subproblem  provides  a basis 
for  making  feasibility  checks  and  could  be  used  to  devise  an  efficient,  implicit 
enumeration  solution  strategy.  Several  approaches  (constraint  aggregation  and 


Generalized  Lagrange  Multipliers)  which  might  be  applied  to  simplify  the  net- 
work design  subproblem  to  a straightforward  flow-circulation  problem  were  in- 
vestigated but  do  not  appear  to  offer  general  computational  advantages.  An 
alternate  model  by  which  tours  could  be  constructed  using  linear  programming 
routines  in  a branch-and-bound  algorithm  provides  a viable  approach.  Contin- 
ued study  and  evaluation  of  this  line  of  investigation  appears  to  offer  the  poten- 
tial for  determining  mathematically  optimal  solutions  to  planning  problems  of 
considerable  size. 

Recommendations  for  Future  Research 

The  most  pressing  needs  at  this  time  are  to  computerize  the  algorithms 
developed  during  this  study  and  to  evaluate  their  effectiveness  in  solving  actual 
problems.  A great  deal  of  work  in  this  direction  has  already  been  accomplished, 
but  the  limitations  imposed  by  the  current  minigrant  have  precluded  complete 
implementation  of  the  computer  program  under  development.  However,  it  is 
believed  that  the  objectives  of  the  minigrant  were  fully  accomplished.  Expecta- 
tions of  computational  efficiency  (based  on  the  experience  with  concepts  proven 
successful  in  many  other  applications)  appear  sufficient  to  conclude  that  the 
approaches  developed  for  the  transportation  system  problem  will  also  be  found 
satisfactory  runtimewise.  An  outline  of  the  existing  computer  program  and 
its  capabilities  is  given  in  Appendix  G. 

In  most  practical  applications  it  is  iinfwrtant,  not  only  to  provide  an 
optimal  solution,  but  also  to  be  able  to  answer  a variety  of  'what  if'  questions  in 
a sensitivity  analysis  designed  to  promote  understanding  of  system  operation. 
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The  approaches  developed  during  this  study  are  readily  amenable  to  such  exten- 


sions.  A full  range  of  such  capabilities  should  be  included  in  any  program  used 
in  an  actual  design  situation. 

Capability  to  formulate  the  complex  systems  planning  problem  and  apply 
efficient,,  large  scale  programming  techniques  should  stimulate  study  of  a range 
of  related  and  extended  problems.  For  example,  the  related  problem  which  re- 
quires minimization  of  (weighted)  ship  time  subject  to  a total  system  cost  re- 
striction involves  a minor  change  in  model  formulation  and  implementation  of 
a modest  number  of  changes  in  the  solution  procedure.  This  problem  is  also 
important  to  the  Logistics  Command  which  must  ultimately  provide  the  best 
level  of  service  possible  within  the  funding  level  appropriated  by  Congress. 
Extended  problems  of  this  general  type  might  involve,  for  example,  assessment 
of  the  ALC  stocking-point  location  policy,  or  evaluation  of  the  location  of  vehicle- 
maintenance  facilities. 
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Appendix  A 

A Summary  of  Current  AFLC 
Operating  Procedures 
(from  the  minigrant  proposal) 


The  Air  Force  expends  over  forty  million  dollars  annually 


to  operate  a dedicated  air  transport  network  called  Logair 
and  an  additional  thirty  million  dollars  for  transporation 
on  commercial  modes.  Even  though  most  of  this  expense  is 
incurred  to  provide  premium  (air)  transportation,  there  is 
no  guarantee  that  the  Order  and  Ship  Times  (0  & ST)  experienced 
are  the  shortest  possible.  While  shipping  time  may  not  consti- 
tute the  bulk  of  0 & ST,  it  does  represent  the  portion  which 
is  most  difficult  to  optimize  from  a systems  viewpoint. 
Transportation  has  a significant  impact  on  the  logistics  system 
not  only  through  direct  costs  but  also  from  the  inventory  invest- 
ment and  operational  readiness  ramifications  of  the  service 
level  provided. 

Material  is  transported  amongst  CONUS  installations  using 
commerical  modes  (such  as  UPS,  Air  Freight,  and  truck)  and  the 
extensive  air  transport  network,  Logair,  which  is  contracted 
especially  for  Air  Force  needs.  Decisions  to  ship  via  a com- 
mercial mode  are  made  on  an  ad  hoc  basis  within  general  guide- 
lines provided  by  the  Air  Force  policy.  A shipment  is  eligible 
for  premium  transportation  on  Logair  if  criteria  specified  by 
management  policy  are  satisfied. 

The  inventory  system,  which  gives  rise  to  transportation 
requirements,  is  composed  of  two  echelons:  the  first  level  is 
distributed  amongst  user  bases  and  the  second  level  for  each 
item  is  located  at  a single  point.  Five  Air  Logistics  Centers 
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(ALC)  in  the  CONUS  provide  world-wide  support  at  this  second 
echelon.  Each  item  is  stocked  at  only  one  ALC.  Furthermore, 
depot  level  repair  for  an  item  is  located  at  the  same  ALC. 
Several  general  types  of  shipments  result  from  this  inventory 
structure . 

Investment  items  probably  constitute  the  bulk  of  ship- 
ments which  require  premium  transportation  on  Logair.  By 
definition,  these  items  are  so  costly  that  repair  is  economi- 
cally justified  in  lieu  of  replacement  upon  failure.  An 
investment  item  may  be  repaired  at  the  local  base  or,  if  nec- 
essary, returned  to  the  depot  at  the  appropriate  ALC  for 
maintenance.  In  the  latter  case,  the  base  would  requisition 
a replacement  from  the  ALC  so  that  the  base  mission  would  not 
be  impared  by  the  reduction  in  the  number  of  items  on  hand. 
This  procedure  corresponds  to  an  (s-l,s)  inventory  policy. 

On  the  average,  some  number  of  items  will  be  in  the  shipping 
pipeline  from  base  to  ALC  and  vice  versa.  Reductions  in 
shipping  time  therefore  have  a proportional  influence  on  the 
number  of  spares  required  to  support  a weapons  system.  Since 
investment  items  are  expensive,  premium  transportation  is 
justified  to  provide  expedited  shipping  times  in  order  to 
reduce  total  system  cost. 

Less  expensive  items  are  controlled  by  an  economic  order 
quantity  (EOQ)  inventory  policy  and  are  not  repaired  upon 
failure.  Base  stock  is  replenished  by  placing  an  EOQ  re- 
quisition on  the  appropriate  ALC.  Management  policy  provides 
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sufficient  safety  stock  to  accommodate  longer  shipping  times 
so  that  more  economical  ground  transport  on  commercial  modes 
may  be  used.  However,  random  demand  patterns  and  other  factors 
may  make  it  necessary  to  upgrade  the  urgency  of  a shipment, 
making  it  eligible  for  air  transportation. 

A variety  of  factors  may  precipitate  critical  situations 
in  which  a weapons  system  would  become  inoperative  (NORS 
for  example)  if  highly  expedited  transportation  is  not  pro- 
vided. In  this  sense,  the  transportation  system  is  relied 
upon  to  compensate  for  complications  which  might  arise  from 
other  aspects  of  the  logistics  system.  Lateral  support  capa- 
bility (base  to  base  resupply)  represents  an  example  of  such 
flexibility  afforded  by  Logair. 

Additionally,  special  shipments  are  made  on  Logair  since 
the  dedicated  mode  provides  certain  flexibilities  not  typically 
provided  by  commercial  modes.  Classified  shipments,  explosive 
and  hazardous  material,  and  outsize  cargo  may  be  shipped  on 
Logair  for  economic  and  managerial  reasons. 

In  all  cases,  the  individual  requisitioner  assigns  a trans- 
portation priority  to  the  shipment.  This  priority  directs 
transportation  personnel  to  provide  either  ground  or  air  service 
and  specifies  a standard  time  within  which  delivery  must  be 
made . 

Shipping  time  is  therefore  the  primary  measure  of  perfor- 
mance for  the  transportation  system.  But  costs  must  be  controlled 
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at  reasonable  levels.  In  general,  the  cost  of  shipping  de- 
creases as  the  quantity  shipped  increases.  For  example,  a 
full  truck  load  ( FTL)  may  be  transported  at  a rate  per  ton- 
mile  which  is  significantly  less  than  that  which  would  be 
incurred  by  several  less  than  truck  load  (LTL)  shipments. 
Additionally,  FTL  shipments  may  provide  better  service  time 
since  cargo  is  not  transferred  amongst  vehicles. 

An  optimal  transportation  plan  should  evaluate  economics 
which  might  result  from  accumulating  individual  shipments  at 
certain  points  in  the  network  in  order  to  achieve  overall  cost 
reductions  from  increasing  vehicle  utilization  and  promoting 
use  of  vehicles  of  suitable  capacities.  While  service  level 
may  be  improved  under  certain  conditions,  a time  penalty  for 
diverting  shipments  to  permit  aggregation  would  be  incurred 
in  most  cases.  The  optimal  plan  must  therefore  achieve  a 
balance  between  shipping  time  performance  and  system  cost. 

Vehicles  which  operate  on  a fixed  schedule  in  a dedicated 
mode  incur  total  costs  which  may  be  independent  of  utilization. 
Therefore,  once  a dedicated  mode  is  established,  shipments 
should  be  consolidated  to  achieve  maximum  vehicle  utilization 
and  to  avoid  additional  costs  from  shipping  by  other  modes. 

The  Logair  cost  structure  includes  mileage  charges  and  a fee 
for  each  landing  so  the  cost  is  independent  of  cargo  carried. 
Non-air  eligible  cargo  may  be  carried  at  times  to  improve 
aircraft  utilization  and  avoid  extra  costs  of  shipping  the 
material  on  a commercial  mode. 
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The  Logair  network  is  designed  with  the  objective  of 
providing  the  service  level  advantages  of  a dedicated  mode. 

The  network  is  specified  by  A FLC  personnel  to  accommodate  the 
requirements  of  CONUS  installations.  Requirements  are  fore- 
cast by  each  installation  and  submitted  for  AFLC  use.  The 
current  design  process  does  not  explicitly  consider  shipping 
time  requirements;  only  total  transportation  capacity  needs 
are  evaluated. 

The  network,  as  historically  designed,  consists  of  trunk 
and  feeder  routes.  Efficient,  high  capacity  aircraft  are  used 
on  the  trunk  lines  to  transport  aggregated  shipments  amongst 
the  five  ALC's,  Wright-Patterson  Air  Force  Base,  and  six 
ports  of  debarkation  from  which  shipments  are  transshipped  to 
overseas  installations . Feeder  routes,  which  are  serviced  by 
smaller  aircraft,  connect  ALC's  to  bases.  Current  management 
policy  permits  only  one  Logair  flight  to  service  a base  each 
day.  Each  of  the  bases  in  the  network  may  actually  represent 
the  collective  needs  of  a number  of  installations  in  its 
vicinity,  so  shipments  may  require  trucking  to  ultimate 
destinations . 

AFLC  personnel  designate  the  aircraft  to  be  used  and  the 
route  each  will  service  by  applying  a heuristic,  manual  de- 
sign procedure.  Since  aircraft  require  extensive  routine 
maintenance,  routes  are  carefully  designed  considering  the 
locations  of  contractor  maintenance  facilities.  The  length 
of  each  route  must  consider  aircraft  flight  capability  as 
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well  as  federal  regulations  regarding  crew  working  times. 

Long  routes  are  feasible  but  increase  costs  since  crews  must 
be  changed. 

The  final  aspect  of  system  design,  specifying  arrival/ 
departure  times  for  each  aircraft  at  each  point  serviced,  is 
currently  negotiated  with  the  contractor.  Prime  considera- 
tions in  this  process  are  the  times  at  which  tases  provide 
ground  crews  to  service  Logair  flights  and  flight  crew  needs 
as  represented  by  the  contractor. 

Management  policies  by  which  the  system  operates  appear 
attractive  intuitively,  but  the  cost  implications  which  they 
entail  have  not  been  well  defined  quantitatively.  Current 
heuristic  procedures  do  not  allow  the  trade-offs  between 
shipping  times  and  costs  to  be  made  explicit  for  management 
review.  Succinctly,  a systems  analysis  is  necessary  to 
prescribe  the  optimal  commercial/dedicated  mode  transportation 
system  and  permit  thorough  evaluation  of  service  level  require- 
ments and  management  policies. 
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Appendix  B 


The  Transportation  System  Planning  Model 
Developed  during  the  1976  USAF/ASEE 
Summer  Faculty  Program 

(Excerpted  from  the  minigrant  proposal) 
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Index  on  points  (ALC’s  and  BASES)  in  the  system 
Index  on  points  (ALC's  and  BASES)  in  the  system 
Index  on  vehicles  (aircraft,  trucks,  commercial 
modes)  available  for  use 

Index  on  transport  chains  connecting  material  flow 

origins  and  destinations 

Cost  for  vehicle  k to  service  point  i 
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r 


^mid 


m 

n 


0 = 
P/m  - 

<5ik  “ 
s = 

t = 

I _ 

( i jk  )m 
TJm  = 
uk  = 
Uim  = 


Vi 


m 


Vk  = 


V 


mid 


X. 

1 jkn 


Y(  i jk)./mn 


Set  of  chains  connecting  origin  i to  destination  d 
for  material  m 

Index  on  material  types  (priority,  weight,  cost,  etc 
Number  of  a specific  arc  in  a vehicle  route  (path) 
or  material  chain 
Point  of  origin  for  a vehicle 

Portion  of  material  m (for  a specific  OD  pair) 
flowing  on  chain  ^ 

1 if  vehicle  k serves  point  i and  "0"  otherwise 
Source  or  origin  of  a material  flow  chain 
Termination  point  for  a vehicle  route  (path) 

Time  required  for  material  m to  traverse  arc  (ijk) 
Time  required  for  material  m to  traverse  chain  ^ 

1 if  vehicle  k is  utilized  and  "0"  otherwise 
Urgency  (or  priority)  factor  associated  with 
material  m on  chain 

1 if  chain^m  is  used  and  "0"  otherwise 
Fixed  cost  for  using  vehicle  k 

Number  of  chains  consddered  for  mid  requirements 

1 if  vehicle  k traverses  arc  (ij)  nth  in  its 

route  (path) , "0"  otherwise 

t h 

1 if  arc  (ijk)  is  nl  in  chain  X for  material 
m,  "0"  otherwise 


(DV)* 

(DV)* 

(DV)* 

(DV)* 

(DV)* 

(DV)* 
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m « ■ Hi  ■ MHl 


= Factor  to  represent  compartmentalized  capacity 
of  vehicles 

Tm  = Maximum  standard  shipping  time  permitted  for 
material  m 

T tcj  = Standard  mean  shipping  time  for  the  system 

(DV)*  Implies  a decision  variable  in  the  program. 
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Transportation  System  Model 
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I . The  Dedicated  Mode  Network 
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, vehicle  k leaves  its  point  of 
origin  0 if  used 

uk  ' Vk'  veh;’-cle  ^ terminates  at  point  t 
if  used 

X.  < u,  ; ¥,  , vehicle  k traverses  its  nth  arc 

l jkn  — k kn 


jtkn 


at  most  once 


jkn  < u^  : k<  vehicle  k traverses  arc  ij  at 

most  once  (optional) 


Ex..k 

. ijkn 
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is  connected 
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(time)  limits 
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(time)  limits  for  arcs  on  tour  k 
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II . 


[Linear  constraints  in  u,  giving  vehicle  relationships] 


Multicommodity  Material  Flow 
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nth  link  at  most  once 


. ...  . /)  < v /i  ;¥../?  v/1  , chain  Am  uses  arc 
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( i jk  )Amn 
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^Y(  jik)/m(n+l)  : *Am  e^mid  > j+t  # j+0 
chain  _/m  is  connected 


Y(  i jk)^mn  — ^—Xijkn  ' Vijk,-/m  , flow  link- 


n " ~mid 

vehicle  arc  relationship 
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Appendix  C 


Second  Model  Formulation  Developed  in  Fall,  1976 
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TIIIO  OHIO  STATE  UNIVERSITY 


December  7,  1976 


Dr.  Ismail  N.  Shimi 
Program  Manager 

Air  Force  Office  of  Scientific 
Research 

Mathematical  and  Information 
Science 

Bolling  Air  Force  Base 
Washington,  D.C.  20332 

Dear  Dr.  Shimi: 

I submitted  a proposal  (number  06017-55-00)  entitled 
"Optimization  of  a Transportation  System  Planning  Problem"  to 
your  minigrant  program  early  this  fall.  The  proposed  research 
is  a continuation  of  work  initiated  during  the  1976  USAF/ASEE 
Summer  Faculty  Research  Program  held  at  the  Logistics  Command 
Headquarters  at  Wright-Patterson  Air  Force  Base.  The  proposal 
presented  a mathematical  programming  formulation  of  a complex 
logistics  problem  and  recommended  research  to  continue  model 
development  and  examine  potential  solution  approaches  to  evaluate 
computation  time  required  to  solve  problems  of  practical  size. 

Continued  thought  on  this  problem  has  just  recently  yielded 
a much  improved  formulation  of  the  transportation  system  pro- 
blem: and,  according  to  our  telephone  conversation  of  12-7-76, 

I am  forwarding  the  new  formulation  to  you  for  review.  While 
this  may  not  yet  represent  the  ultimate  form,  it  is  a significant 
improvement  over  the  model  in  the  proposal,  and  I believe  it 
should  be  considered  as  a part  of  the  proposal  review  process. 

The  new  model  contributes  to  the  first  step  of  the  proposed 
work  by  providing  a more  tractable  form  and  by  making  a number 
of  feasibility-checking  procedures  obvious.  The  new  model  makes 
a significant  reduction  in  the  number  of  zero/one  integer  variables 
and  should,  therefore,  reduce  computation  times  required  to  solve 
problems  of  realistic  size.  A portion  of  the  network  design 
burden  is  switched  from  the  master -integer  program  to  new  variables 
Y and  W which  may  be  treated  as  continuous  variables.  In  fact, 
specific  values  of  Y and  W are  either  zero  or  one  and  are  deter- 
mined directly  by  current  values  of  the  integer  variables  X and 
G in  constraints  (7),  (8),  (9),  and  (19).  Values  of  dual  variables 
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Dr.  Ismail  N.  Shimi 
December  7,  1976 
Page  two 

associated  with  these  constraints  may  be  determined  directly 
(this  portion  of  the  problem  is  trivial  and  does  not  even  require 
a linear  programming  solution)  and  provide  feedback  to  improve 
the  network  design  in  the  X and  G variables. 

Constraints  in  the  X variables  (equations  2-5)  represent  a 
special  case  transportation  problem  as  do  those  constraints 
involving  only  the  G variables  (equations  12-15).  These  two  sets 
of  constraints  are  tied  together  by  equation  (18).  This  structure 
permits  efficient  feasibility  checking,  and,  in  fact,  may  lead  to 
an  efficient  solution  algorithm. 

This  formulation  is  unique  and  provides  an  efficient  approach 
to  modeling  a class  of  problems  in  which  a network  must  be  designed 
as  a component  of  the  overall  problem.  Problems  in  scheduling, 
sequencing,  and  vehicle  routing  fall  into  this  class  and  certain 
problems  not  amenable  to  previous  approaches  may  be  solvable  using 
this  decomposition  approach.  A formal  research  proposal  to  study 
the  structure  of  this  class  of  problems  is  being  written  and  will 
be  submitted  to  you  in  the  near  future.  The  transportation  system 
problem  is  one  specific,  albeit  complex,  model  in  this  class. 

I would  appreciate  your  considering  the  improved  formulation 
as  a factor  in  the  review  of  the  minigrant  proposal  submitted 
earlier. 


Sincerely , 


W.  E.  Wilhelm 
Assistant  Professor 


WEW : cm 

cc:  J.  Abell,  AFLC 

A.  B.  Bishop,  OSU 
M.  W.  Woody,  OSURF 
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The  Transportation  Systems  Problem 


Notation 


A'/ 

A.  . . 
ljnkq 

Ck 

c>(F/e > 


dijk 

Fijk 

Dk'-k 


Dk'§k 


~msd 


i 

j 

k 


Administrative  cost  of  material  flow  chain  £ 

Cost  of  material  flow  chain  q £ using  arc  ijnk 

Weight  (volume)  capacity  of  vehicle  k 

Cost  of  shipping  amount  of  material  on  chain  jl 

Destination  point  for  material  shipment 

Cost  of  vehicle  k traversing  arc  ij 

Distance  on  arc  ij  for  vehicle  k 

Upper,  lower  bounds  on  total  route  (or  path) 
distance  (time)  for  vehicle  k 

Upper,  lower  bounds  on  distance  (time)  for 
each  arc  traversed  by  vehicle  k 

Required  weight  (volume)  of  material  m to 
be  shipped  from  S to  d 

Weight  (volume)  flow  of  material  m on  chain  X 

1 if  material  flow  chain  H departs  from  point 
i as  the  qth  link  and  0 otherwise  (DV*) 

Index  on  points  in  the  system 

Index  on  points  in  the  system 

Index  on  vehicles  (aircraft,  trucks,  commercial 
modes,)  available  for  use 


K 


Number  of  vehicles  available  for  use 


/ 


L 


ik 


Index  on  transport  chains  connecting  material 

flow  origins  and  destinations.  Each  k is 

specific  to  some  f , . 

msd 

Cost  for  vehicle  k to  service  point  i 
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dC.  msd 


Set  of  chains  connecting  source  s to  destination 
d for  material  m 


I 


m 


N, 


[Nvl 


0 

P 


Pink  St 


q 

Q/ 

[Q  ^ 3 


Tink/ 

T/ 

!nk 

VJL 

V/q 

V,. 


Vmsd 

X . , 
ink 


Index  on  material  types  (priority,  weight, 
cost , etc . ) 

Number  of  a specific  arc  in  a vehicle  route 
( path ) 

Number  of  points  which  vehicle  k may  visit 
n 1 , 2 , • • • , 

set  of  specific  N^  specific  points  which  k may  visit 

Point  of  origin  for  a vehicle 

Portion  of  f , which  flows  on  chain 
msd 

P^  on  arcs  (ink)  used  by  chain  Jl 

Index  on  number  of  arcs  in  material  flow  chain 

Maximum  number  of  arcs  management  will  permit  in 
chain  / 

Set  of  points  which  might  be  used  in  chain  / 

Large  positive  number 

Source  or  origin  of  a material  flow  chain  Z 
Termination  point  for  a vehicle  route  (path) 

Time  required  to  traverse  arch  (ink)?) 

Time  required  for  material  m to  traverse  chain  . fl- 
1 if  vehicle  services  n points,  0 otherwise 
Urgency  (or  priority)  factor  associated  with  chain  / 
1 if  chain  / is  composed  of  q arcs,  0 otherwise 
Fixed  cost  for  using  vehicle  k 

Number  of  chains  considered  for  msd  requirements 

til 

1 if  point  i is  the  n point  in  the  route  (path) 
for  vehicle  k 0 otherwise 
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(DV)  * 
(DV)* 


(DV)* 

(DV)* 

(DV) 


1 if  vehicle  k travels  from  point  i to  point  j 
0 otherwise 


(DV)* 


1 if  both  G and  X allow  vehicle  k and  chain  ji 
to  traverse  the  arc  from  point  i to  point 
j,  0 otherwise 

Factor  to  represent  compartmentalized  capacity 
of  vehicles 

Maximum  standard  shipping  time  permitted  for 
material  m 

Standard  mean  shipping  time  for  the  system 


(DV)*  Implies  a decision  variable  in  the  program. 


Yijk 


Wi  jnkqJ" 


m 


std 
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(DV)* 


min 


[c1  (P.f  ,)  + Z d'  • Y . . 

. £ ' je  msd  • i]K  ID 


i jk 


jk  i jk 


+ I hk  <£Xink>  Vk  (1-U0k»  + le  A'(  U-V0^  > 

ik  n k * 


+ Z Ai  jnkqj?.  Wi  jnkq^ 
lkl 


Subject  to: 


I.  Dedicated  Mode  Network  Constraints 


? Xink  - (E  n Unk>'  Vk 
in  n 


Z nx.  v = 7Z  (n2+n)  U , , Vk 


Link  2 
m n 


xWk  = ( 1_U0k)  ' Vk 


nX 


= (Z  nU  , ) for  k which  traverse  a path; 


tnk  " vn  “ nk 


n=2 , 3 , . . . <N, 


N, 


Z Unk 


U , = 1 , Vk 


n=0 


X01k  + Xj2k  - Y0jk  ^ 1 Vfc  je[V  ' W 


Xink  + Xj,n+l,k  Yijk 


1 Vk;ije[Nk],  i^j 
n=2 , 3 , . . . ( N^-l 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


-nX  . , + y'nU  . + Y > 1 for  k which  traverse 

:nk  L nk  3?,n+l  - a route. 


je[N,],  j/0 : n=2, 3,  . . . ,N, 


-k  - £ Dijk  Yijk  - Dk  , 


Bk^Dijk  Yijk^gk  , Vk‘  ^^Nk] 


II.  Multicommodity  Material  Flow  Chains 


% Giq a = ^ qvq  je  v/ 

lq  M q ^ 


£ ^iqt"  2 q Vq pVt 


GslJt~  (1_v0j?)  ' 


q Gdq^  = ^ q Vqj2^  ' q 


E vq/g=  1 . V/ 


l 

^CoCmsd 


(1"V0^  } ^ Vmsd  ' V msd 
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Giqi’+  Gj,q+l,j?  - ^ ^Xink  + Xj,n+l,k^'  Vijq/ 


o + G • . i o + x-  i + X . , _ W.  . . . < 

iq X ],q+l,k  ink  j,n+l,k  ljnkq^  — 


Wijnkq£—  Giq/  ' Wi  jnkq,^  — G j ,q+l , J?  ' Wijnkq^  — Xink 


l jnkq  / — X j , n+1 , k 


Vijql  nk 


III.  Multicommodity  Material  Flow  Constraints 


msd  f 


— 1.0  , Vmsd 


£ lW&  Gsl£ 


6 ink*  «•  E wijnkq^  .V  kn,i* 


Z 6ink*  ^ 1 ■ 


Pink<>£  *T  ° ^-‘inkj?  > Vi 


1 
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IV 


Service  Level  Constraints 


1 


I f , P.  , „ < a,  C,  . V ink 
^ msd  inki  — k k' 

Z U T.  , / f , P.  . n < T , f , ,V  msd  (26) 

. , / ink  H msd  inkX  — msd  msd 

ink 

Z E/msd 

Tink  « fmsd  Pinkj?  - Tstd  ^ fmsd  (27) 

ink/  X msd 

V.  Non-negativity  Constraint 

X, G,6 ,U,v=  0,1  (28) 

Y, W,P,  P > 0 


(25) 


A 


The  objective  of  this  mathematical  programming  formu- 
lation (Equation  1)  is  to  minimize  total  system  cost  which 
is  composed  of  the  cost  of  shipping  on  each  material  flow 
chain,  the  travel  cost  of  vehicles  used  by  dedicated  modes, 
the  cost  of  a vehicle  to  service  a point,  a fixed  charge 
for  each  vehicle  used,  an  administrative  fee  for  each 
material  flow  chain  established,  and  the  cost  of  material 
flow  on  each  arc  in  the  transportation  network. 

The  purpose  of  each  constraint  is  described  below: 

Constraint  Equation 

Number Affect  of  Constraint 


2 assures  that  each  vehicle  services 
the  assigned  number  of  points,  (nU^). 

3 assures  that  vehicle  tours  are  connected 


4 


5 

6 

7 


8 

9 


assures  that  a vehicle  tour  begins  at 
the  predefined  point  of  origin  if  the 
vehicle  is  used 

requires  a vehicle  path  to  terminate 
at  the  predefined  termination  point 

determines  the  number  of  points  in 
the  tour  for  each  vehicle 

defines  the  arc  which  a vehicle  takes 
upon  departing  its  predefined  point  of 
origin 

defines  the  arcs  used  on  the  tour  of 
vehicle  k 

defines  the  arc  from  the  (nUn^)t^1 

point  visited  to  the  predefined  point 
of  origin  for  vehicles  which  traverse 
a route 
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Constraint  Equation 

Number Affect  of  Constraint 


10  limits  vehicle  route/path  distance/time 
to  predetermined  specification 

11  limits  vehicle  arc  distance/time  to 
predetermined  specification 

12  assures  that  each  material  flow  chain 
passes  thru  the  assigned  number  of 
points,  (q Vo) 

Hi 

13  assures  that  material  flow  chains  are 
connected 

14  requires  a material  flow  chain  to 
originate  at  the  source  point 

15  requires  a material  flow  chain  to  end 
at  the  destination  point 

16  determines  the  number  of  points  in  a 
material  flow  chain 

17  assures  that  the  actual  number  of  material 
flow  chains  defined  for  an  (mid)  combina- 
tion is  within  a limit  set  by  management 

18  assures  that  a material  flow  link 
from  point  i to  point  j is  not  defined 
unless  some  vehicle  services  that  arc 

19  defines  the  specific  vehicles  (k) 

which  are  available  for  use  by  the  material 
flow  link  (ijqf  ) from  point  i to  point  j 

20  assures  that  all  material  for  a given 
(mid)  combination  will  be  delivered  to 
the  destination 


21 

22 


allows  material  to  flow  only  in  chains 
which  are  defined 

the  vehicle  used  by  chain  X to  depart 
point  i must  be  one  which  traverses 
the  correct  arc  (note:  this  constraint 
holds  for  fixed  values  of  W and  carries 
implications  about  off  loading  material 
on  vehicles  which  visit  a point  or 
an  arc  more  than  once. ) 
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Constraint  Equation 
Number 


Affect  of  Constraint 


23 

24 


25 

26 


27 


a material  flow  link  will  use  at  most 
one  vehicle  between  points  i and  j 

assures  that  all  vehicle  arcs  used  by 
material  flow  chain  £ carry  the  same 
material 

assures  that  vehicle  capacity  is  not 
violated  on  any  vehicle  arc 

assures  that  an  urgency/material 
weighted  time  average  of  shipping  time 
for  each  (mid)  combination  in  within 
limits  set  by  management 

assures  that  an  urgency/material  weighted 
time  average  of  shipping  time  over  all 
material  flow  chains  is  within  a standard 
set  by  management  for  overall  system 
performance 
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Appendix  D 

Summary  of  Notation  Used  in  Report 
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I 


A 


b 

q 


R 

i 

B-1 


B 


O 


r rA 
Lmsd 

crT 

msd 


C 


rk 


cl 


d 


set  of  extreme  points  in  flow  circulation  network  defined  in 
Problem  MP 


1 if  vehicle  k traverses  arc  q (from  point  i to  point  j),  0 otherwise 

1 if  flow  chain  JL  uses  arc  q (ij  and  k implicit),  0 otherwise 

(vehicle)  capacity  provided  on  arc  q 
q =>  specific  points  ij  and  vehicle  k 

The  index  q also  indicates  the  1st  or  2nd  time  arc  ij  is  traversed  by  k. 


transshipment  capacity  at  point  i 

inverse  of  matrix  whose  columns  correspond  to  a basic  feasible 
solution 

minimum  cost  of  making  all  material  shipments  on  commercial  modes 
subject  to  the  ship  time  constraints 

non-negative  constant  in  Benders'  cut  ic  (see  equation  [38]  ) 
set  of  arcs  beginning  at  point  i for  vehicle  k 
cost  coefficient  for  flow  on  chain  SL 
fixed  cost  for  using  vehicle  k 

cost  of  traversing  arc  (ijk)  and  of  servicing  point  j 
cost/unit  to  ship  material  combination  msd  on  commercial  air 
cost/unit  to  ship  material  combination  msd  on  commercial  truck 
cost  of  tour  r for  vehicle  k 

^ Crk  Trk’  Trk  £ \ 
destination  point  for  a material  shipment 

dual  variable  associated  with  constraint  defining  P^  in  figure  2 
length  in  distance  (or  time)  of  arc  q 
maximum  tour  length  for  vehicle  k 
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D 


= dummy  destination  point 


r 


= set  of  tours  excluded  at  step  1/B  in  solving  Problem  01 
D(N)  = 'distance'  from  source  S to  node  N 
Dq  = 'distance'  or  length  of  arc  q 

- set  of  subscripts  of  tours  which  have  no  associated  y values 
Ej(k)  = set  of  arcs  ending  at  point  i for  vehicle  k 

fmsd  = s^iPPinS  requirements  commodity  m,  from  source  s to  destination  d 
F = Problem  PI  objective  function  value 

<r  . = 1 if  flow  chain  H accommodates  combination  msd  = n,  0 otherwise 

i = 1,  2, . . . , I points 

j = 1,  2, ....  I points 

J = vector  of  the  T , variables  that  have  been  assigned  a binary 

specification  inme  current  partial  solution  in  Problem  01 

J - number  of  feasible  vehicle  network  designs 

k = 1,  2, . . . , K vehicles 

k - vehicle  which  is  the  complement  of  vehicle  k 

= column  representing  flow  chain  Ji 
K ' - number  of  vehicle  tours  in  a complete  network 

JL  = 1,  2, . . . , L flow  chains 

m - 1,  2, . . . , M commodities  types  (M  may  imply  all  types) 

n|v.  = number  of  tours  defined  in  Problem  01  for  vehicle  k 

Ntc  - set  of  subscripts  of  the  T variables  which  are  candidates  for 
improving  the  current  partial  solution  in  Problem  01 

O order  and  ship  time  at  source  s for  material  type  m 

sm 

O 'dummy'  point  of  origin  for  all  vehicle  tours 

6 point  of  origin  for  vehicle  k 
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q 


r 


^msd 


q 

s 

m 


Si 

o rA 
msd 

SrT 

smsd 


t. 


i l 


rp  e 

li. 


T-S 

i,  m 


ok 


T 

T 


T 


T(k) 

Tq(k) 


portion  of  material  of  type  t,  fmscj*  shipped  on  flow  chain  JL 
!i  and  t=>  a specific  msd  combination 

Problem  LP  objective  function  value 

1,  2, ....  Q vehicle  arcs  (ijk  implicit) 

index  for  vehicle  tours  in  Problem  PI 

artificial  variable  for  material  flow  requirement  f 

source  point  for  a material  shipment 

dummy  source  point 

slack  variable  associated  with  vehicle-arc  constraint  q 

slack  variable  associated  with  ship  time  constraint  m 

slack  variable  associated  with  transshipment  point  i capacity 

portion  of  f ^ shipped  on  commercial  air 

portion  of  fmS([  shipped  on  commercial  truck 

1 if  flow  chain  SL  transships  at  point  i,  0 otherwise 

shipment  time  on  vehicle  arc  q 

(overnight)  ship  time  delay  incurred  if  flow  chain  X is  carried 
through  vehicle  k termination  point  i 

transshipment  time  for  material  type  m at  point  i 

1 if  vehicle  k is  not  used,  0 if  it  is  used  (the  null  tour)  (T^  - 1 - TQ|.) 

'dummy'  termination  point  for  all  vehicle  tours 

termination  point  for  vehicle  k 

vector  of  Tj,^  values 

set  of  tours  which  vehicle  k might  traverse 
set  of  tours  for  vehicle  k which  traverse  arc  q 
1 il  vehicle  k traverses  tour  r,  0 otherwise 
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u 


number  of  arcs  on  a tour 


= unload  time  incurred  at  shipment  destination  d for  material  type  m 
= set  of  all  vehicles  provided 

= set  of  subscripts  of  vehicles  which  have  been  assigned  (one 
Trk  = 1 in  Jt>  111  Jt 

= set  of  subscripts  of  vehicles  which  are  excluded  because  their  com- 
plementary vehicles  are  in  the  set 

= set  of  vehicles  which  are  candidates  for  inclusion  in  solution 
at  level  one 

set  of  vehicles  for  which  all  tours  C Ut 
= set  of  vehicles  for  which  all  tours  C 

set  of  vehicles  which  have  all  tours  at  level  zero  in  J, 


Crk  - Z b I y I the  coefficient  of  tour  r for  vehicle  k 
QeTk(q)  1 in  Benders'  cut  k 


Z W . „ T . for  T , C J 
kr  rk/c  rk  rk  - t 


1 if  flow  chain  JL  is  carried  through  vehicle  k termination  point  i, 
0 otherwise 

amount  of  'flow'  through  point  i in  subnetwork  k in  Problem  P2 

= continuous  variable  in  Benders'  subproblems 

current  best  Z value  resulting  from  a complete  solution 

= Z value  which  results  from  a completion  of  the  current  partial 
solution,  J. 


parameter  0 £ a . — 1 


dual  variable  associated  with(material  m)ship  time  constraint 


a 1 o 
m ' 1 ivl 


dual  variable  associated  with  transshipment  capacity  at  point  i 


X 

X 

X 


p (m)  = 

^o 

^x 


>7 


T 


T 


T 


(t)  _ 

j i 

A 

msd 

T 

msd 


dual  variable  associated  with  vehicle  arc  q 

ship  time  standard  for  material  type  m 

dual  variable  associated  with  n = msd  ship  requirement 

see  equation  [18] 

iteration  number  between  Benders'  subproblems 

current  Benders'  iteration  number 

continuous  decision  variable,  0 - V - 1 

set  of  flow  chains  accommodating  material  type  m 

dual  variable  associated  with  normalizing  constraint 

dual  variable  associated  with  Benders'  cut  x 

dual  variable  associated  with  tour  length  constraint  for  vehicle  k 

set  of  vehicles  used  simultaneously  in  a shortest  path  problem 
used  for  column  generation 

ship  time  on  flow  chain  JL  (t  implicit) 

ship  time  for  msd  combination  on  commercial  air 

ship  time  for  msd  combination  on  commercial  truck 
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Matrix  Notation  and  Dimensions 


Column 

Reference 

Matrix 

Vector 

Elements 

Dimension(s) 

Figure 

A 

0,1:  to  relate  points  to  a tour  (#  points,  # tours) 

14 

Af 

*msd  Aq£ 

(L,  Q) 

2 

B 

Benders' cut  constants,  B 

K 

/c 

14 

crA 

c rA 

msd 

msd 

2 

CrT 

crT 

'“msd 

msd 

2 

crA 

rA 

msd  msd 

msd 

2 

c rT 

f c rT 

msd  msd 

msd 

o 

[C-by] 

tour  cost  in  Benders'  cut,  ecj. 

[43]  # tours 

14 

G 

gn,£ 

(msd,  L) 

2 

I 

identity  matrix 

appropriate 

2,14 

£ 

right  hand  side  in  network 

logic  constraints, 

eq.  [42] 

14 

p 

p/> 

L 

2 

rb 

artificial  variables,  eq.  [43] 

14 

rk 

artificial  variables,  eq.  [78] 

K 

14 

Rr 

artificial  variables,  Rmscj, 

msd 

2 

eq?S[12] 

sa 

sa 

q 

Q 

2 

SB 

surplus  variables,  eq. |43[ 

K. 

14 

S«£ 

slack  variables,  eq.  [42] 

appropriate 

14 

sp 

slack  variables,  eq.  [40] 

K 

14 

srA 

srA, 

msd 

msd 

2 

SrT 

srT 

“msd 

msd 

2 
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1 

Column 

Reference 

Matrix 

Vector 

Elements 

Dimension  (s) 

Figure 

ss 

Ss 

M 

2 

sl 

s‘ 

1 

# transshipment  points 

2 

T 

Trk 

# tours 

14 

r 

'm 

M 

14 

A 

0, 1:  relate  tours  to  network 

(#  network  constraints. 

14 

constraints 

# tours) 

nt 

0, 1:  relate  flow  chains  to 

(#  transshipment 

2 

transshipment  points 

points,  L) 

Z 

0, 1:  relate  tours  to  vehicles 

(K,  # tours) 

14 

Tf 

fmsd  TV 

(M,  L) 

2 

A 

Tf 

f A 

Imsd  T msd 

(M,  msd) 

2 

T 

Tf 

f t T 

msd  msd 

(M,  msd) 

2 

i 

column  vector  of  ones 

appropriate 

2,14 
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Appendix  F 

Proof  of  an  Assertion  Used  in  Chapter  IV 


Assertion: 


y =0  for  arcs  q not  used  by  the  vehicle  network  specified 
at  iteration  ic  in  LP. 

Proof:  Suppose  the  variables  are  renumbered  so  that  the  first  M are  in  solution. 
T he  row  would  be  of  the  form 


Aj  P - 0 

and  the  first  M columns  would  contain  zeroes,  since  no  ilow  chain,  X , 
which  uses  arc  q may  be  used  in  a feasible  solution.  Currently, 


B (P  SrA  SrT)'=  b,  (P  SrA  SrT)'=  B-1b,  Pq  CR(P  SrA  Sr  1 ), 

(a  P y d0  S ) = CR  B_1  A. 

Add  row  q to  the  problem  and  consider  the  possibility  of  adding  an  addi- 
tional basic  variable,  P . 


and  dual  variables  are  defined  by 
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^ - — — - JMILJLJ 


cbb_1  \ = / S®'1 

CPq  / ' 0 

since  P is  a slack  with  C = 0 
q Pq 

The  dual  variable  associated  with  arc  q is  therefore  zero  as  originally  stated. 
This  argument  may  easily  be  extended  to  prove  the  assertion  for  any  number 
of  constraints  not  made  explicit  in  LP. 
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Appendix  G 
Computer  Program 
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A computer  program  to  test  the  solution  procedures  described  in  Chapters 
III  and  IV  was  written  and  is  currently  being  debugged.  The  program  includes  a 
Revised  Simplex  algorithm  with  Generalized  Upper  Bounding  which  has  been 
specialized  to  improve  run  time  and  reduce  storage  requirements.  The  column 
generation  routine  defined  in  Chapter  III  is  included.  Techniques  such  as  'checks 
for  zeroes'  and  explicit  storage  of  nonzero  elements  only  are  incorporated  as  is 
the  use  of  double  precision  to  moderate  accumulation  of  numerical  round-off 
errors.  Routines  which  implement  the  tour  selection  procedure  (the  algorithm 
for  Problem  01)  have  also  been  keypunched,  along  with  the  necessary  input,  out- 
put, and  cross-communication  routines.  Programs  to  solve  the  two  subproblems 
consist  of  3233  cards  (several  input/ output  routines  were  duplicated  to  allow 
parallel  debugging  of  the  program  halves). 

The  objective  of  promoting  run  time  efficiency  guided  program  design 
throughout.  Other  features  were  included  to  simplify  implementation.  Examples 
are  the  use  of  one-dimensional  storage  vectors  to  reduce  run  time  and  to  allow 
the  routines  to  be  stored  on  disc.  Dimensions  of  arrays  need  be  defined  only  in 
the  main  program  for  each  different  application.  In  addition,  the  program  was 
designed  to  allow  easy  testing,  including  provisions  for  inputting  17  and  values  a 
(see  Chapters  III  and  IV,  respectively)  as  data  elements. 

Emphasis  was  placed  on  pursuing  the  line  of  investigation  reported  in 
Chapter  V in  lieu  of  concentrating  on  the  computer  program.  Investigation  of 
tour  construction  ability  was  deemed  more  important,  overall,  in  this  develop- 
ment study.  Additionally,  this  lint?  of  investigation  was  emphasized  considering 


the  related  research  proposed  as  another  project  and  under  review  by  the  USA  F/ 
OSR  during  most  of  this  summer. 

Variable  definitions  and  program  codes  are  included  in  this  appendix  for 
completeness.  Faculty  release  time  provided  this  fall  quarter  which  is  not  used 
in  report  preparation  will  be  applied  to  continue  computer  program  development. 
However,  the  level  of  testing  required  should  be  accomplished  on  a project  of 
longer  duration  and  more  concentrated  effort. 
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LIST  OF  VARIABLES 


COLSUM R)  ( 1 Wr| .^column  sum  in  Benders’  cuts  NTOUR 


o 

o 


Variable  Definition  Int  * 2 Real*8  M ax  # rows  M ax  # col 


NMATST  number  of  (material  categories)  ship  time  constraints 

NMAT  + 1 =>  all 


Variable  Definition  Int  * 2 Real'S  Max#  rows  Max#  col 


r 
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NTH  variable  number  of  P - MSD  + MTH 


Variable  : Definition  Int  * 2 Real *8  Max#  rows  Max#  col 


CODES  FOR  SUBROUTINE 


AIESGE  (II,  12,  13, 14, 15,  16, 17,  18, 19,  110,  111,  112,  113,  114,  115) 


11  CURRENT  LP  SOLUTION 

KEYVAR(.  ),  NET,  P(.  ) 

BASIC (.  ),  KNET(.  ),  P(.  ) 

K VAR(.  ) 

12  BINV(.  ) 

13  ALP(.  ) , LALP(.  ) 

14  ENTERING  GENERATED  COLUMN  COL(.  ) 

15  CUAIN(.  ) 

16  AS(.  ),  DBARS(.  ),  ABARS(.  ) 

17  DUAL(.  ) 

18  W(.  ),  WMIN(. ) 

19  01  SOLUTION  STATUS 

T(.  ),  JT(.  ),  NOEK(.  ),  V(.  ),  LB,  ZC,  UB 

110  NODE(.  ) 

11 1 CUMCST(.  ),  PATH(.  ),  ORDERl(.  ),  ORDER2(.  ) 

112  Reserved  tor  Problem  MP 

113 

114  " 

115  CURRENT  01  DECISION  CRITERIA 

WSUM(.  ),  CTMIN(. ),  COLDEC(.  ),  COLMAX(.  ),COLSUM(.  ) 
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A 


ERROR  CODES 


PROGRAM 

SUBROUTINE 

SUBROUTINE  NAME 

2 

1 

READIN 

2 

2 

PREPLP 

2 

3 

RSGUB 

2 

4 

COLGEN 

2 

5 

SHORT 

2 

6 

REORDR 

2 

7 

PHILP 

2 

9 

UPBINV 

3 

1 

READIN 

3 

2 

PREP01 

3 

3 

IE01 

3 

4 

PREPND 

PROGRAM  CODES 


2 subroutines  associated  with  problem  LP 

3 subroutines  associated  with  problem  01 
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A 


ROW  AND  VARIABLE  INDICES  FOR  PROBLEM  LP 

ROW  INDICES  NUMBER 


(1) 

Ship  time  constraints 

NMATST 

(2) 

Transshipment  capacity  constraints 

N TRANS 

(3) 

Vehicle  arc  capacity  constraints 

QMAX 

(4) 

PQ  row  defining  the  LP  objective 

1,  the  MTH  row 

VARIABLE  INDICES 

NUMBER 

(1) 

Artificial  variables  for  ship  requirements 

MSD 

(2) 

Slacks  for  ship  time  constraints 

NMATST 

(3) 

Slacks  for  transshipment  capacity  constraints 

N TRANS 

(4) 

Slacks  for  vehicle  arc  capacity  constraints 

QMAX 

(5) 

Pq  the  LP  objective  function  value 

1,  the  NTH  VARIABLE 

(6) 

Commercial  air  shipment  for  msd  requirement 

MSD 

(7) 

Commercial  truck  shipment  for  msd  requirement 

MSD 

(8) 

Generated  columns 

up  to  MSD  + MTH  - 1 
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